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Abstract 


This  report  describes  a  comprehensive  study  of  the  statistical  characteristics  of  concentra¬ 
tion  fluctuations  in  a  neutrally  buoyant  tracer  plume  dispersing  through  a  large  array  of 
building-like  obstacles,  each  of  which  measured  12.2  m  x  2.42  m  x  2.54  m.  The  plumes 
were  released  both  upwind  and  within  the  obstacle  array  for  a  range  of  source  heights  be¬ 
tween  0.15  and  5.2  m.  Detailed  flow  field  and  instantaneous  plume  concentration  data  were 
obtained  from  a  comprehensive  series  of  tracer  experiments  which  utilized  a  large  number 
of  high-resolution  concentration  detectors,  accompanied  by  the  simultaneous  acquisition  of 
meteorological  and  turbulence  measurements  with  sonic  anemometer/thermometers.  Ex¬ 
tensive  analyses  are  performed  on  the  plume  concentration  data,  and  results  are  presented 
for  a  number  of  concentration  statistics  such  as  the  mean  plume  lateral  and  vertical  spreads, 
mean  concentration,  fluctuation  intensity,  peak  concentration  to  concentration  standard  de¬ 
viation  ratio,  concentration  probability  density  function  (pdf),  concentration  power  spectra, 
and  various  concentration  time  and  length  scales  of  dominant  motions  in  the  array  plume 
(e.g.,  integral  scale,  Taylor  microscale). 

For  the  range  of  downwind  distances  from  the  source  examined,  the  lateral  mean  concentra¬ 
tion  profiles  are  well  approximated  by  a  Gaussian  distribution.  The  vertical  profiles  of  mean 
concentration  develop  in  a  rather  complex  manner  with  downwind  distance,  with  the  result 
that  the  reflected  Gaussian  form  is  generally  a  less  than  ideal  description  of  the  mean  array 
plume  in  the  vertical  direction.  A  comparison  of  the  array  plume  with  an  open-terrain 
plume  as  a  function  of  downwind  distance  indicates  that  the  obstacle  array  significantly 
increases  the  lateral  and  vertical  plume  spreads  and  decreases  the  magnitude  of  the  plume 
centreline  mean  concentration.  The  small-scale,  high-intensity  turbulence  generated  in  the 
obstacle  array  results  in  a  drastic  reduction  in  the  concentration  fluctuation  level  in  the 
array  plume  compared  to  an  open-terrain  plume  under  similar  conditions.  The  evolution 
of  the  concentration  pdf  at  a  fixed  range,  but  with  decreasing  height  from  above  and  into 
the  obstacle  array  is  similar  to  that  obtained  at  a  fixed  height  but  with  increasing  down¬ 
wind  distance  from  the  source.  The  integral  and  Taylor  microscale  time  and  length  scales 
of  the  plume  increase  significantly  within  the  obstacle  array.  Concentration  power  spectra 
measured  within  the  array  had  a  greater  proportion  of  the  total  concentration  variance  in 
the  lower  frequencies  (energetic  subrange),  with  a  correspondingly  smaller  proportion  in 
the  higher  frequencies  (inertial-convective  subrange).  It  is  believed  that  these  effects  re¬ 
sult  from  the  rapid  and  efficient  stirring  and  mixing  of  plume  material  by  the  small-scale, 
high-intensity  turbulence  within  the  array. 


DRDC  Sufficld  TR  2003-097 


1 


Resume 


Ce  rapport  decrit  uiie  etude  approfondie  des  caracteristiques  statistiques  des  fluctuations  de 
la  concentration,  dans  un  panache  traceur  d’une  densite  de  flottabilite  neutre,  se  dispeisant 
a  travel's  un  important  reseau  d'obstacles  ressemblants  a  des  batiments,  chacun  inesmant 
12,2  m  X  2,42m  x  2,54  m.  Les  panaches  ont  ete  disperses  centre  le  vent  et  a  I  inteiieur 
du  reseau  d’obstacles,  a  une  hauteur  variant  de  0.15  a  5,2  m.  Des  champs  de  propagation 
detailles  et  des  donnees  instantanees  de  concentrations  de  panaches  ont  ete  obtenus  a  partir 
de  series  completes  d’essais  a  partir  de  traceurs  utilisant  un  grand  nombre  de  detecteurs 
de  concentration  de  haute  resolution.  Ceci  tout  en  effectuant  I’acquisition  simultanee  de 
mesures  meteorologiques  et  des  turbulences,  au  moyen  d’anemometres  et  thermometres 
soniques.  Des  analyses  extensives  ont  ete  effectuees  sur  les  donnes  de  concenti ation  du 
panache  et  les  resultats  sont  pr&entes  pour  un  certain  nombre  de  statistiques  de  concenti  a- 
tion  telles  que  les  ecarts  lateraux  et  verticaux  du  panache  moyen,  la  concentration  moyenne, 
I’intensite  de  la  fluctuation,  le  rapport  concentration  maximale  —  ecart-type  de  la  concen¬ 
tration  maximale,  la  fonction  de  densite  de  la  concentration,  les  spectres  de  la  puissance  de 
la  concentration  ainsi  ejue  les  echelles  des  durees  et  longueurs  variees  de  concentration  des 
motions  dominantes,  a  I’interieur  du  panache  du  reseau  (par  ex  :  1  echelle  integrale.  1  echelle 
microscopique  de  Taylor). 

En  ce  ejui  concerne  la  portee  des  distances,  sous  le  vent  a  partir  de  la  source,  qui  ont  ete 
examinees,  les  profils  lateraux  de  la  concentration  moyenne  sont  bien  estimfe  par  une  loi  de 
Gauss.  Les  profils  verticaux  de  concentration  moyenne  se  developpent  d  une  manieie  assez 
complexe.  dans  la  distance  sous  le  vent,  le  resultat  etant  que  la  forme  gaussienne  ciui  est 
reflechie  forme  une  description  du  panache  du  reseau,  dans  le  vertical,  qui  n  est  pas  ideale. 
Une  comparaison  du  panache  du  reseau  avec  ceux  ejui  se  developpent  dans  un  champ  libie  en 
fonction  de  la  distance  sous  le  vent,  indique  ejue  le  r&eau  d’obstacles  augmente  de  manieie 
signifiante  I’etendue  du  panache  lateral  et  vertical  et  diminue  I’ampleur  de  la  concentration 
de  la  ligne  mediane  du  panache.  La  turbulence,  a  petite  echelle  et  de  haute  intensite.  generee 
dans  le  reseau  d’obstacles,  resulte  en  une  reduction  drastique  du  niveau  de  fluctuation  de  la 
concentration  dans  le  panache  du  rseau,  compare  celui  en  champ  libre,  dans  des  conditions 
similaires.  L’evolution  de  la  fonction  de  densite  de  la  concentration  a  une  distance  donnee. 
mais  d’une  hauteur  decroissante  a  partir  du  dessus  et  vers  I’interieur  du  reseau  d  obstacles, 
est  similaire  a  celle  obtenue  a  une  hauteur  definie  mais  a  une  distance  sous  le  vent  croissante 
a  partir  de  la  source.  Les  echelles  de  duree  et  celles  de  longueur  integrals  et  a  echelles 
microscopiques  de  Taylor  du  panache  augmentent  de  maniere  signifiante  a  I’interieur  du 
reseau  d'obstacles.  Les  spectres  de  puissance  de  la  concentration  mesures  a  l  interieur  du 
r&eau  avaient  une  proportion  plus  importante  de  la  variance  de  la  concentration  totale 


n 


DRDC  Sufficld  TR  2003-097 


dans  les  frequences  plus  faibles  (sous  intervalles  d’energie)  avec  une  proportion  plus  petite 
correspondante  dans  les  plus  hautes  frequences  (sous  intervalles  par  inertie  -  de  convection). 
On  pense  que  ces  effets  resultent  du  brassage  et  du  melange  des  matieres  du  panache  causes 
par  la  turbulence,  a  faible  echelle  et  de  haute  intensite,  a  rinterieur  du  r&eau. 


hi 


DRDC  Suffield  TR  2003^097 


Executive  summary 


Introduction:  It  is  anticipated  that  Canadian  Forces  (CF)  in  the  foreseeable  future  will 
have  to  fight  in  or  protect  urban  areas,  whether  in  battle,  peace-making,  peacekeeping, 
or  counter-terrorist  operations.  The  inci'eased  awareness  and  importance  accorded  b>  the 
public  worldwide  and  their  governments  to  maintain  appropriate  defences  against  chemical 
and  biological  warfare  (CBW)  agents  in  an  urban  (built-up)  environment,  the  prediction 
of  casualties  and  human  performance  degradation  resulting  from  such  releases,  and  the 
development  of  operational  procedures  and  regulations  to  control,  mitigate,  and  monitor  the 
fate  of  CBW  agents  in  urban  areas  with  high  population  densities,  will  require  mathematical 
modeling  of  urban  wind  flows  and  dispersion.  One  of  the  drawbacks  in  the  computer 
simulation  of  urban  flow  and  dispersion  is  the  lack  of  experimental  data.  This  capability 
gap  impairs  the  development  and  thorough  validation  of  numerical  models.  That  difficulty 
has  been  partly  addressed  recently  in  the  case  of  the  design  and  implementation  of  the  Mock 
Urban  Setting  Trial,  for  which  extensive  field  measurements  have  been  made  available  for 
the  development  and  validation  of  urban  flow  and  dispersion  models. 

Results:  This  report  describes  a  comprehensive  study  of  the  statistical  characteristics  of 
concentration  fluctuations  in  a  neutrally  buoyant  tracer  plume  dispersing  through  a  large  ar¬ 
ray  of  building-like  obstacles,  each  of  which  measured  12.2  m  x  2.42  m  x  2.54  m.  The  plumes 
were  released  both  upwind  and  within  the  obstacle  array  for  a  range  of  source  heights  be¬ 
tween  0.15  and  5.2  ni.  Detailed  flow  field  and  instantaneous  plume  concentration  data  were 
obtained  from  a  comprehensive  series  of  tracer  experiments  which  utilized  a  large  numbei 
of  high-resolution  concentration  detectors,  accompanied  by  the  simultaneous  acquisition  of 
meteorological  and  turbulence  measurements  with  sonic  anemometer/thermometers.  Ex¬ 
tensive  analj'ses  are  performed  on  the  plume  concentration  data,  and  results  are  presented 
for  a  number  of  concentration  statistics  such  as  the  mean  plume  lateral  and  vertical  spreads, 
mean  concentration,  fluctuation  intensity,  peak  concentration  to  concentration  standard  de¬ 
viation  ratio,  concentration  probability  density  function  (pdf),  concentration  power  spectra, 
and  various  concentration  time  and  length  scales  of  dominant  motions  in  the  airay  plume 
(e.g.,  integral  scale.  Taylor  microscale). 

For  the  range  of  downwind  distances  from  the  source  examined,  the  lateral  mean  concentra¬ 
tion  profiles  are  well  approximated  by  a  Gaussian  distribution.  The  vertical  profiles  of  mean 
concentration  develop  in  a  rather  complex  manner  with  downwind  distance,  with  the  result 
that  the  reflected  Gaussian  form  is  generally  a  less  than  ideal  description  of  the  mean  array 
plume  in  the  vertical  direction.  A  comparison  of  the  array  plume  with  an  open-terrain 
plume  as  a  function  of  downwind  distance  indicates  that  the  obstacle  array  significantly 
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increases  the  lateral  and  vertical  plume  spreads  and  decreases  the  magnitude  of  the  plume 
centreline  mean  concentration.  The  small-scale,  high-intensity  turbulence  generated  in  the 
obstacle  array  results  in  a  drastic  reduction  in  the  concentration  fluctuation  level  in  the 
array  plume  compared  to  an  open-terrain  plume  under  similar  conditions.  The  evolution 
of  the  concentration  pdf  at  a  fixed  range,  but  with  decreasing  height  from  above  and  into 
the  obstacle  array  is  similar  to  that  obtained  at  a  fixed  height  but  with  increasing  down¬ 
wind  distance  from  the  source.  The  integral  and  Taylor  microscale  time  and  length  scales 
of  the  plume  increase  significantly  within  the  obstacle  array.  Concentration  power  spectra 
measured  within  the  array  had  a  greater  proportion  of  the  total  concentration  variance  in 
the  lower  frequencies  (energetic  subrange),  with  a  correspondingly  smaller  proportion  in 
the  higher  frequencies  (inertial-convective  subrange).  It  is  believed  that  these  effects  re¬ 
sult  from  the  rapid  and  efficient  stirring  and  mixing  of  plume  material  by  the  small-scale, 
high-intensity  turbulence  within  the  array. 

Significance  and  Future  Plans:  The  primary  purpose  of  this  report  is  to  present  a 
detailed  picture  of  the  concentration  statistics  in  plumes  from  point  sources  dispersing 
through  a  large  array  of  obstacles.  The  aim  is  to  provide  valuable  basic  data  on  the  struc  ture 
and  development  of  array  plumes  that  can  be  used  to  guide  the  construction  and  assessment 
of  predictive  models  for  dispersion  in  a  built-up  area. 


Yee,  E.  2003.  Mock  Urban  Setting  Trial:  Data  Analysis  and  Interpretation.  DRDC  Suffield 
TR  2003-097,  Defence  R&D  Canada  -  Suffield. 
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Sommaire 


Introduction:  On  prevcit  que.  dans  un  avenir  assez  rapproche,  les  Forces  canadiennes 
(CF)  devront  combattre  a  Tinterieur  des  zones  urbaines  ou  bien  les  defendre,  soit  an  cours 
de  batailles  soit  an  cours  d'operations  de  maintien  de  la  paix  ou  anti-terronstes.  L’opmion 
publique  et  les  gouvernements  du  monde  entier,  mieux  avertis,  accordent  de  plus  en  plus 
d’importance  a  la  necessity  de  maintenir  un  niveau  approprie  de  defense,  centre  les  agents 
de  guerre  chimiques  et  biologiques  (CBW),  dans  les  milieux  urbains  (batis).  On  prevoit 
aussi  que  le  nombre  de  victimes.  la  degradation  des  performances  humames  resultant  de 
telles  emissions  ainsi  que  la  mise  au  point  des  procedures  operationnelles  et  des  reglements 
d’exploitation  mis  en  place  pour  controler,  contenir  et  surveiller  le  devenir  des  agents  de 
guerre  CB  dans  les  zones  urbaines,  ayant  une  population  de  haute  densite,  necessitera  une 
modelisation  mathematique  du  mouvement  des  vents  et  celle  de  la  dispersion  d’un  polluant 
atmospherique.  Un  des  inconvenients  de  la  simulation  par  ordinateur  des  mouvements  et  de 
la  dispersion  en  zone  urbaine  est  le  manque  de  donnees  experimentales.  Cette  lacune  com- 
promet  la  mise  au  point  et  la  validation  complete  des  modeles  num^-iques.  Cette  difficulte 
a  ete  partiellement  resolue  recemment  dans  le  cas  de  la  conception  et  de  I’lmplementation 
d’un  Essai  fictif  en  milieu  urbain  pour  lequel  ont  ete  prises  des  mesures  elaborees  sur  le 
terrain  afin  de  realiser  la  mise  au  point  et  la  validation  des  modeles  de  mouvements  et  de 
dispersion  en  milieu  urbain. 

Resultats:  Ce  rapport  decrit  une  etude  approfondie  des  caract&istiques  statistiques  des 
fluctuations  de  la  concentration,  dans  un  panache  traceur  d’une  densite  de  flottabflite  neu- 
tre,  se  dispersant  a  travers  un  important  reseau  d’obstacles  ressemblants  a  des  batiments, 
chacun  mesurant  12, 2  m  x  2. 42m  x  2, 54  m.  Les  panaches  ont  ete  disperses  centre  le  vent 
et  a  I’interieur  du  reseau  d'obstacles,  a  une  hauteur  variant  de  0,15  a  5,2  m.  Des  champs 
de  propagation  detailles  et  des  donnees  instantanees  de  concentrations  de  panaches  ont  ete 
obtenus  a  partir  de  series  completes  d’essais  a  partir  de  traceurs  utilisant  un  grand  nombre 
de  detecteurs  de  concentration  de  haute  resolution.  Ceci  tout  en  effectuant  I’acquisition 
simultanee  de  mesures  meteorologiques  et  des  turbulences,  au  moyen  d’anemometres  et 
thermometres  soniques.  Des  analyses  extensives  ont  t  effectues  sur  les  donnes  de  concen¬ 
tration  du  panache  et  les  rsultats  sont  prsents  pour  un  certain  nombre  de  statistiques  de 
concentration  telles  que  les  ecarts  lateraux  et  verticaux  du  panache  moyen,  la  concentration 
moyenne,  I’intensite  de  la  fluctuation,  le  rapport  concentration  maximale  —  ecart-type  de 
la  concentration  maximale,  la  fonction  de  densite  de  la  concentration,  les  spectres  de  la 
puissance  de  la  concentration  ainsi  que  les  echelles  des  durees  et  longueurs  variees  de  con¬ 
centration  des  motions  dominantes,  a  I’interieur  du  panache  du  reseau  (par  ex  :  I’echelle 
integrale,  I’echelle  microscopique  de  Taylor). 
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En  ce  qui  concerne  la  portee  des  distances,  sous  le  v^ent  a  partir  de  la  source,  qui  out  ete 
examinees,  les  profils  lateraux  de  la  concentration  moyenne  sont  bien  estimes  par  une  loi  de 
Gauss.  Les  profils  verticaux  de  concentration  moyenne  se  developpent  dtme  maniere  assez 
complexe.  dans  la  distance  sous  le  vent,  le  resultat  etant  que  la  forme  gaussienne  qui  est 
reflechie  forme  une  description  du  panache  du  reseau,  dans  le  vertical,  qui  n’est  pas  ideale. 
Une  comparaison  du  panache  du  reseau  avec  ceux  cjui  se  developpent  dans  un  champ  libre  en 
fonction  de  la  distance  sous  le  vent,  indique  que  le  reseau  d’obstacles  augmente  de  maniere 
signifiante  letendue  du  panache  lateral  et  vertical  et  diminue  I'ampleur  de  la  concentration 
de  la  ligne  mediane  du  panache.  La  turbulence,  petite  echelle  et  de  haute  intensite,  generee 
dans  le  reseau  d’obstacles,  r&ulte  en  une  reduction  drastique  du  niveau  de  fluctuation  de  la 
concentration  dans  le  panache  du  rseau,  compare  celui  en  champ  libre.  dans  des  conditions 
similaires.  L'evolution  de  la  fonction  de  densite  de  la  concentration  a  une  distance  donnee. 
mais  d’une  hauteur  decroissante  a  partir  du  dessus  et  vers  I’interieur  du  reseau  d’obstacles. 
est  similaire  a  celle  obtenue  a  une  hauteur  definie  mais  a  une  distance  sous  le  vent  croissante 
a  partir  de  la  source.  Les  echelles  de  duree  et  celles  de  longueur  integrales  et  a  echelles 
microscopiques  de  Taylor  du  panache  augmentent  de  maniere  signifiante  a  I’interieur  du 
reseau  d'obstacles.  Les  spectres  de  puissance  de  la  concentration  mesures  a  I’interieur  du 
reseau  avaient  une  proportion  plus  importante  de  la  variance  de  la  concentration  totaJe 
dans  les  frequences  plus  faibles  (sous  intervalles  denergie)  avec  une  proportion  plus  petite 
correspondante  dans  les  plus  hautes  frequences  (sous  intervalles  par  inertie  -  de  convection). 
On  pense  que  ces  effets  resultent  du  brassage  et  du  melange  des  matieres  du  panache  causes 
par  la  turbulence,  a  faible  echelle  et  de  haute  intensite.  a  I’interieur  du  reseau. 

La  portee  des  resultats  et  les  plans  futurs:  Le  but  principal  de  ce  rapport  est  de 
presenter  en  details  les  statistiques  de  concentration  dans  les  panaches,  a  partir  de  leur 
point  d'origine  de  dispersion  et  a  travers  un  important  reseau  d’obstacles.  L’objectif  est  de 
fournir  des  donnes  fondamentales  precieuses  au  sujet  de  la  structure  et  du  developpement 
du  panache  d’un  reseau  qui  pen  vent  etre  utilisees  comme  guide  dans  la  construction  et 
revaluation  de  modeles  de  prdiction  de  la  dispersion  dans  des  zones  baties. 
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Introduction 


In  recent  years,  owing  to  an  increasing  likelihood  of  an  accidental  or  deliberate  release  of 
a  hazardous  material  in  an  urban  (built-up)  area  such  as  an  industrial  complex  or  residen¬ 
tial  and  commercial  dense  urban  estate  where  the  population  density  is  high,  considerable 
attention  has  focused  on  the  development  of  models  to  predict  the  dispersion  of  toxic  mate¬ 
rials  in  the  urban  environment.  Many  studies  of  urban  dispersion  in  the  past  have  restricted 
attention  to  the  far-held  of  the  toxic  release  where  the  plume  dimensions  are  signihcantly 
larger  than  the  individual  buildings.  In  this  case,  the  enhanced  dispersion  arising  from 
the  (usually  complicated)  conhguration  of  buildings  in  the  urban  complex  is  characterized 
simply  as  an  increased  surface  roughness  leading  to  a  greater  intensity  of  turbulence  for 
plume  dispersion  (Martin  [1];  Gifford  and  Hanna  [2]).  However,  it  is  expected  that  the 
highest  concentrations  and  the  greatest  potential  for  deleterious  effects  of  exposure  to  the 
toxic  material  will  occur  in  the  near-held  of  the  toxic  release.  In  consequence,  a  number 
of  researchers  (e.g.,  Meroney  [3]  :  Hosker  [4]:  Hosker  and  Pendergrass  [5])  have  contributed 
to  the  important  topic  of  the  how  and  near-held  dispersion  around  a  single  building  or  a 
small  group  of  buildings.  Nevertheless,  owing  to  its  intrinsic  complexity,  the  problem  of 
near-held  plume  dispersion  through  large  groups  of  obstacles  has  received  attention  only  in 
recent  years.  In  particular,  the  systematic  investigations  initiated  by  Baechlin  et  al.  ([6], 
[7]),  Jerram  et  al.  [8],  Davidson  et  al.  ([9],  [10]),  and  Macdonald  et  al.  ([11],  [12])  are 
notable  and  have  contributed  signihcantly  to  our  current  understanding  of  the  behaviour 
of  a  passive  neutrally-buoyant  plume  passing  through  a  large  array  of  obstacles.  Hanna 
and  Chang  [13]  describe  the  Kit  Fox  roughness  obstacle  held  experiments  conducted  at  the 
Fi-enchman  Flat  area  of  the  Nevada  Test  Site  in  late  August  and  early  September,  1995. 
These  experiments  were  designed  to  represent  an  oil  rehnery  or  chemical  plant  at  approxi¬ 
mately  1:10  scale,  and  included  some  ground-level  dense  gas  releases  with  the  initial  cloud 
depth  less  than  the  obstacle  height.  In  this  study,  the  authors  noted  that  the  peak  mean 
concentration  in  the  cloud  increased  with  increasing  cloud  advection  speed,  and  observed 
that  this  relationship  appeared  to  vary  with  the  obstacle  packing  density  which  in  turn 
dictated  whether  the  increased  turbulence  intensity  or  decreased  mean  wind  speed  in  the 
obstacle  array  produced  the  dominant  effect. 

Virtually  all  the  research  effort  on  the  dispersion  of  plumes  through  groups  of  obstacles  has 
focused  exclusively  on  the  description  of  the  mean  concentration  field  (the  concentration 
field  is  described  only  in  a  time-averaged  sense).  The  only  notable  exceptions  to  this  are 
afforded  by  the  study  of  Davidson  et  al.  [9]  who  investigated  the  effect  of  an  obstacle  array 
on  the  structure  of  the  plume  from  a  continuous  point  .source  and  by  the  investigation  of 
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Pavag6au  and  Schatzmann  [14]  who  exantiined  the  statistical  propeities  of  the  instantaneous 
concentration  field  produced  by  a  continuous  line  source  within  a  street  canyon  flow.  The 
paucity  of  information  concerning  the  statistical  description  of  the  concentration  fluctua¬ 
tions  of  a  plume  as  it  passes  through  a  large  group  of  obstacles  is  a  rather  surprising  state 
of  affairs.  Indeed,  concentration  fluctuations  are  a  ubiquitous  feature  of  dispersing  plumes 
and,  in  recent  years,  it  has  been  recognized  that  this  feature  has  piactical  importance  in  a 
number  of  engineering  and  technological  applications.  Some  examples  include  the  proba¬ 
bility  of  ignition  of  inflammable  gases  which  depends  upon  the  instantaneous  concentration 
lying  between  the  lower  and  upper  flammability  limits,  irrespective  of  the  mean  concentra¬ 
tion;  the  a,ssessment  of  human  health  hazard  from  toxic  substances  where  the  response  is 
nonlinear;  and,  estimation  of  the  perception  of  odours  which  depends  on  whether  instanta¬ 
neous  concentrations  are  likely  to  exceed  certain  critical  threshold  levels  over  short  periods 
of  time. 

The  importance  of  concentration  fluctuations  coupled  with  the  recent  availability  of  fast- 
response  instrumentation  for  the  measurement  of  the  phenomenon  has  led  to  various  exper¬ 
imental  studies  that  have  contributed  significantly  to  the  extensive  body  of  knowledge  that 
now  exists  concerning  the  statistical  properties  of  concentration  fluctuations  in  dispersing 
plumes.  These  studies  have  been  mostly  limited  to  dispersion  over  flat  and  relatively  smooth 
terrain  in  order  to  minimize  the  complexities  of  the  analysis  and  have  included  measure¬ 
ments  made  in  wind  tunnels  and  water  channels  (e.g..  Fackrell  and  Robins  [15];  Stapountzis 
et  al.  [16];  Bara  et  al.  [17]  as  well  as  full-scale  atmospheric  measurements  (e.g.,  Sawford  et 
al.  [18];  Dinar  et  al.  [19];  Peterson  et  al.  [20];  Mylne  and  Mason  [21];  Mylne  [22];  Yee  et  al. 
[23];  Yee  et  al.  [24],  [25];  among  others). 

Because  full-scale  atmospheric  studies  of  the  detailed  structure  of  fluctuating  concentrations 
in  a  plume  dispei'sing  through  a  large  array  of  building-like  obstacles  are  rare,  we  have 
conducted  a  comprehensive  series  of  tracer  experiments  with  the  goal  of  providing  more 
extensive  measurements  of  concentration  fluctuation  statistics  in  an  array  plume.  The 
Mock  Urban  Setting  Trial  (MUST)  described  herein  was  designed  to  provide  insight  on  the 
instantaneous  dispersion  of  a  tracer  through  a  large  array  of  building-like  obstacles.  The 
present  study  w’as  motivated,  in  part,  by  a  need  to  provide  a  better  understanding  of  how 
the  structure  of  a  plume  is  modified  as  it  disperses  through  a  large  array  of  obstacles  and 
to  subsequently  identify  the  physical  mechanisms  responsible  for  this  modification.  With 
this  in  mind,  mean  concentration  field  characteristics  and  higher-order  moment  properties 
of  the  instantaneous  plume  concentration  are  measured  at  a  large  number  of  positions 
downwind  of  point  sources  located  both  upwind  of  and  within  an  obstacle  array  foi  a 
number  of  different  release  heights.  Even  though  the  mean  concentration  field  will  be 
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studied,  this  paper  is  concerned  primarily  with  the  effects  of  a  large  group  of  obstacles  on 
plume  concentration  fluctuations,  as  it  is  in  this  area  that  the  main  virtues  of  the  work  lie. 
It  is  believed  that  the  present  field  experiment  (viz.,  MUST)  is  a  logical  addition  to  the 
series  of  previous  experimental  works  on  mean-field  dispersion  in  array  plumes  and  on  the 
study  of  concentration  fluctuations  in  open-terrain  plumes. 
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Site,  topography,  and  obstacle  array 


Site  and  topography 


The  measurements  were  carried  out  in  September  2001  at  Horizontal  Grid  on  U.S.  Army 
Dugway  Proving  Ground  which  is  located  in  the  Great  Basin  Desert  of  northwestern  Utah, 
40°  12.606'  N,  113°  10.635'  W.  The  site  elevation  is  1310  m  above  mean  sea  level.  The 
test  site  is  predominantly  flat.  The  cross-sections  in  contour  maps  of  the  area  surrounding 
the  test  site  show  that  the  terrain  slopes  gently  upwards  to  the  south  with  a  slope  of 
approximately  0.0005.  Terrain  features  that  may  influence  winds  over  Horizontal  Grid 
include  sand  dunes  4  to  6  m  in  height  located  at  about  1  km  to  the  north.  At  about  12  km 
southeast  of  the  site.  Granite  Mountain  rises  700  m  above  the  basin  floor.  Approximately 
24  km  northeast  of  the  site  is  Cedar  Mountain  whose  ridge  rises  600  m  above  the  basin 
floor. 

During  the  experiments,  the  horizontally  homogeneous  site  was  covered  with  a  mixture  of 
sparse  greasewood  and  sagebrush  ranging  in  height  from  about  0.4  to  0.75  m.  The  average 
aerodynamic  roughness  length,  zq,  and  the  displacement  height,  d,  which  were  determined 
from  mean  wind  profiles  measured  under  near-neutral  stratification  (where  the  mean  wind 
speed  variation  with  height  can  be  represented  by  a  simplified  semi-logarithmic  relation) 
were  4.5  ±0.5  cm  and  0.37 ±0.09  m,  respectively.  Both  zq  and  d  were  not  dependent  on  wind 
direction.  The  diagnosed  value  of  d  is  somewhat  uncertain,  but  appears  consistent  with  the 
common  rule  of  thumb  that  the  displacement  height  should  be  chosen  to  be  approximately 
three-quarters  of  the  height  of  the  roughness  elements.  The  diagnosed  value  of  zq  appears 
to  be  relatively  insensitive  to  the  exact  value  of  d. 


Obstacle  array 


Each  obstacle  was  a  rectangular  parallelepiped,  with  a  width  {W)  of  12.2  m.  length  (L) 
of  2.42  m,  and  height  (H)  of  2.54  m.  A  total  of  120  obstacles  was  placed  in  an  aligned 
configuration  consisting  of  12  rows  of  10  obstacles  each  with  an  average  obstacle  spacing  of 
(u')  =  7.9  m  in  the  spanwise  (lateral,  or  y-)  direction  and  (/)  =  12.9  m  in  the  streamwise 
(or.  X-)  direction  (see  Figure  1).  Consequently,  the  overall  width  and  length  of  the  obstacle 
array  were  193  m  and  171  m,  respectively,  giving  a  spanwise  aspect  ratio  of  width-to-height 
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{Ly/H)  of  75.9  and  a  lengthwise  aspect  ratio  of  length-to-height  {Lx/H)  of  67.3  for  the 
array.  This  gives  a  frontal  area  density  (frontal  area  of  obstacle  (W  x  H)  divided  by  the  lot 
area  per  obstacle)  of  0.10  and  a  plan  area  density  (plan  area  of  obstacle  {L  x  IT)  divided  by 
the  lot  area  per  obstacle)  of  0.096  for  the  obstacle  array.  For  ease  of  reference,  each  row  and 
column  of  the  obstacle  array  was  assigned  a  letter  from  A  through  L  and  a  number  from  0 
through  9,  respectively  (Figure  1).  In  this  convention,  the  obstacle  at  the  intersection  of  row 
A'  and  column  k  will  be  denoted  Xk  where  X  €  {^4,  .B,  C, . . . ,  L}  and  fc  G  {0, 1. 2, . . . ,  9}. 


The  normal  to  the  long  face  of  the  obstacles  in  the  array  was  oriented  30°  west  of  north 
to  take  advantage  of  the  prevailing  wind  directions  at  the  test  site  in  the  sector  from  SSW 
to  SE,  resulting  in  a  flat  and  homogeneous  upwind  fetch  of  more  than  10  km.  The  x 
(streamwise,  or  longitudinal)  and  y  (spanwise,  or  transverse)  coordinates  are  defined  so 
that  a:  =  0  is  the  front  edge  of  the  obstacle  array  and  y  =  0  is  the  midpoint  (centreline) 
of  the  obstacle  array.  The  z  (vertical)  coordinate  is  defined  so  that  2  =  0  is  the  ground 
surface  beneath  the  obstacles.  Consequently,  with  this  coordinate  system,  the  positive 
.c-axis  is  directed  30°  west  of  north,  whereas  the  positive  y-axis  is  directed  30°  south  of 
west.  The  unit  vectors  along  the  x-,  y-.  and  2-axes  form  a  right-handed  coordinate  frame. 
Furthermore,  the  following  notation  is  used  throughout  this  study:  ui  =  u.  U2  =  v,  uz  =  w 
are  the  instantaneous  velocity  components  along  the  x-,  y-,  and  2-directions,  respectively; 
and,  Ui  and  u[  are  the  mean  and  turbulent  fluctuations  of  the  velocity,  respectively,  with 
[/,.  =  Ui  -b  u[  (hence,  u\  is  the  deviation  of  the  total  velocity  Uj  from  its  mean  \alue  Ui).  In 
this  paper,  a  primed  variable  will  be  used  to  denote  a  departure  of  that  variable  from  its 
mean  value. 
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Instrumentation  and  experimental  details 


Details  of  the  instrumentation  deployed  and  the  experiments  conducted  in  MUST  are  given 
in  Biltoft  [26].  so  only  a  brief  summary  of  the  instrumentation  used  to  acquire  the  data 
analysed  here  vill  be  presented. 


Mean  velocity  and  turbulence  measurements 


Reliable  measurements  of  the  undisturbed,  upwind  flow  that  can  be  used  to  determine 
the  state  of  the  equilibrium  surface-layer  flow  approaching  the  obstacle  array  is  essential. 
Measurements  of  the  vertical  profiles  of  the  mean  horizontal  wind  velocity  and  turbulence 
in  the  upwind  flow  were  obtained  from  a  16-m  telescoping  pneumatic  mast  erected  30.5  m 
(12JT)  upwind  of  the  midpoint  (centreline)  of  the  front  of  the  obstacle  array.  This  mast  was 
instrumented  with  three  horizontal  2-dimensional  (2-D)  sonic  anemometer/thermometers 
(Applied  Technologies  Inc.,  Model  RSWS-211/2SX,  15  cm  path  length)  at  the  4-,  8-,  and 
16-m  levels.  Similarly,  2-D  sonic  anemometer/thermometers  (Applied  Technologies  Inc.. 
Model  RS'V\’S-211/2SX,  15  cm  path  length)  were  mounted  at  the  4-,  8-,  and  16-m  levels 
of  a  16-m  pneumatic  mast  deployed  30.5  m  {\2H)  downwind  of  the  back  of  the  obstacle 
array.  The  sonic  anemometer/thermometers  were  oriented  30°  west  of  true  north  to  take 
advantage  of  the  prevailing  winds  within  the  SSW  to  SE  sector. 

Measurements  of  mean  velocity  and  turbulence  were  obtained  also  within  and  above  the 
obstacle  array.  Vertical  profiles  of  mean  wind  speed  and  turbulence  were  obtained  from 
the  32-m  lattice  tower  located  near  the  centre  of  the  obstacle  array  (Figure  1).  This  tower 
was  instrumented  with  four  3-dimensional  (3-D)  sonic  anemometer/thermometers  (Applied 
Technologies  Inc.,  Model  RSWS-211/3SX,  15  cm  path  length)  at  the  4-,  8-,  16-,  and  32- 
m  levels.  Furthermore,  four  6-m  lattice-type  towers  were  erected  at  or  near  the  centre 
of  the  four  quadrants  of  the  array.  These  four  towers  (Figure  1),  which  were  labeled  as 
towers  A  (NW  quadrant),  B  (SW  quadrant),  C  (NE  quadrant),  and  D  (SE  quadrant)  were 
each  instrumented  with  3-axis  Solent  WindMaster  ultrasonic  anemometers  (Gill  Instruments 
Ltd)  at  the  2.4-  and  6-m  levels.  Finally,  four  3-D  sonic  anemometer/thermometers  (Applied 
Technologies  Inc.,  Model  SATI/3VX,  15  cm  path  length)  were  deployed  at  a  height  of  1.15  m 
within  a  street  canyon  to  study  the  pattern  of  the  flow  between  two  rows  of  obstacles. 
Three  of  these  sonic  anemometers  were  equally  spaced  in  the  streamwise  direction  across 
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the  centre  of  the  street  canyon  between  obstacles  and  G6  (Figure  1),  whereas  the 
fourth  sonic  anemometer  was  centred  in  the  spanwise  gap  between  obstacles  G6  and  G7. 
It  should  be  pointed  out  that  the  probe  (VX;*  of  these  sonic  anemometers  was  designed 
primarily  for  turbulence  measurements  in  canopies  where  the  wind  speed  is  low  and  the 
wind  direction  highly  unpredictable.  The  Applied  Technologies  Inc.  and  Gill  Windmaster 
sonic  anemometer  data  were  recorded  at  10  Hz  and  20  Hz,  respectively. 


Concentration  measurements 


In  each  experiment,  propylene  (CsHe)  tracer  gas  was  released  at  a  point  from  a  specially 
designed  gas  dissemination  system  that  used  a  mass  flow  controller  to  maintain  a  constant 
flow  rate.  The  gas  dissemination  system  consisted  of  three  propylene  cylinders  connected 
in  parallel  and  immersed  in  a  hot-water  bath.  A  Matheson  lL-510  regulator  was  connected 
to  the  outlet  of  the  cylinders  to  ensure  a  constant  downstream  pressure,  and  a  flexible  hose 
was  used  to  connect  the  regulator  to  the  inlet  fitting  of  the  mass  flow  controller  (Teledyne 
Hastings- Raydist).  This  controller,  which  consisted  of  a  sensor,  electronic  circuitry,  a  shunt, 
and  a  valve,  was  used  to  set,  control,  and  measure  the  flow  rate  of  gas  through  the  dis¬ 
semination  system.  The  controller  allowed  a  constant  gas  flow  rate  to  be  maintained  at 
a  user-selected  reference  level  between  20  and  225  1  min  ^  to  within  about  2%.  A  quick- 
release  connector  mated  the  outlet  fitting  of  the  mass  flow  controller  to  the  dissemination 
hose,  which  was  connected  to  the  base  of  the  disseminator,  a  schedule  40  PVC  (polyvinyl 
chloride)  pipe  1  m  in  length  and  0.05  m  in  diameter  through  which  the  gas  was  released. 

The  instantaneous  plume  concentrations  were  measured  with  two  types  of  fast-response 
photo-ionization  detectors;  namely,  digital  photo-ionization  (dPID)  detectors  manufactuied 
commercially  by  Aurora  Scientific  Inc.  (Aurora.  Ontario,  Canada)  and  Ultra  Violet  Ion 
Collector  (UVIC)  detectors  made  by  Industrial  Development  Bangor  Ltd  (United  Kingdom). 
These  detectors  draw  a  stream  of  air  into  a  sampling  volume  that  is  exposed  to  a  stream 
of  ultra-i  iolet  (UV)  radiation  with  a  photon  energi'  of  about  10.6  eV,  which  is  sufficient  to 
ionize  the  tracer  gas  propylene.  The  ions  produced  are  collected  on  an  electrode  system,  the 
current  being  converted  to  a  voltage  signal  which  is  digitized  using  a  PC-controlled,  high¬ 
speed,  high-resolution  16-bit  analog-to-digital  (A-D)  input-output  board  and  multiplexer. 
Both  types  of  detectors  provide  a  frequency  response  of  50  Hz  with  a  sensitivity  of  about 
0.01  ppm  by  volume  of  propylene. 

Horizontal  profiles  of  concentration  statistics  were  measured  using  40  dPIDs  which  were 
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arrayed  along  four  horizontal  sampling  lines  that  were  parallel  to  and  centred  in  the  street 
canyons  (Figure  1).  Sampling  lines  1.  2.  3.  and  4  were  centred  in  the  street  canyons  between 
obstacle  rows  I  and  J.  G  and  H,  E  and  F,  and  C  and  D.  respectively;  and,  these  sampling 
lines  consisted  of  12,  9.  9,  and  10  detectors,  respectively.  The  concentration  detectors  along 
the  four  horizontal  sampling  lines  were  placed  at  a  height,  2^.  of  1.8  m  (2^/^  =  0.71). 
Vertical  profiles  of  concentration  statistics  were  characterized  by  8  dPIDs  deployed  on  the 
32-m  lattice  tower  near  the  centre  of  the  obstacle  array  at  heights  of  1,  2,  4,  6,  8,  10,  12, 
and  16  m  and  by  24  UVICs  mounted  on  the  four  6-m  towers  A,  B,  C,  and  D.  On  each  of 
these  6-m  towers,  6  UVICs  were  deploj-ed  at  the  following  levels:  1,  2,  3,  4,  5,  and  5.9  m. 

In  each  experiment,  the  propylene  gas  was  released  and  sampled  continuously  over  a  period 
of  approximately  15  min.  At  the  flow  rates  for  the  gas  release  (between  150  and  225  1  min~^), 
the  negative  buoyancy  effects  of  the  gas  were  insignificant  (the  release  was  passive).  Time 
series  of  instantaneous  concentrations  were  measured  at  various  points  in  both  horizontal 
and  vertical  cross-sections  through  the  dispersing  plume  at  various  downwind  distances 
from  the  source,  and  for  five  source  heights  (2*  =  0.15,  1.3.  1.8.  2.6,  and  5.2  m;  or,  z^/H  = 
0.059,  0.51,  0.71,  1.02,  and  2.05).  Measurements  were  made  for  sources  positioned  upwind 
of,  within,  and  above  the  obstacle  array.  A  total  of  63  continuous  source  tracer  release 
experiments  were  attempted.  Twenty-one  selected  experiments  were  analyzed  in  detail. 
Pertinent  tracer  release  parameters  from  these  selected  field  experiments  are  summarized  in 
Table  1.  In  Table  1,  Q  is  the  tracer  release  rate  at  the  source,  and  (x^,  j/j,  Zg)  is  the  location 
of  the  source  with  respect  to  the  coordinate  system  defined  in  Figure  1.  The  start  time  for 
each  experiment  refers  to  the  time  the  gas  tracer  was  first  released  from  the  source. 


Data  reduction 


The  velocity  time  series  from  the  sonic  anemometer/thermometers  were  examined  for  out¬ 
liers  and  missing  data  points.  Values  of  the  instantaneous  velocity  that  exceeded  the  limits 
of  four  standard  deviations  from  the  mean  level  were  removed  from  the  time  series.  The 
removed  and  missing  data  points  were  replaced  using  linear  interpolation.  However,  out¬ 
liers  and  missing  data  points  in  the  velocity  time  series  were  rare  (occurring  less  than 
about  0.0001%  of  the  time).  The  corrected  velocity  time  series  from  the  sonic  anemome¬ 
ter/thermometers  were  rotated  in  the  horizontal  plane  to  bring  the  u-  and  u-components  of 
the  velocity  vector  in  alignment  with  horizontal  axes  aligned  along  the  x-  and  y-directions 
(Figure  1)  of  the  obstacle  array.  Linear  trends  were  removed  from  these  time  series  before 
calculation  of  relevant  variances  and  covariances.  This  allowed,  at  each  3-D  sonic  anemome- 


DRDC  Sufficld  TR  2003-097 


9 


tei'/thermometer  probe,  the  nieasurement  of  eddy  correlation  stresses  and  heat  fluxes,  as 
well  as  the  mean  wind  speed,  mean  wind  direction,  and  temperature.  The  relative  impor¬ 
tance  of  shear  and  buovancy  to  the  turbulence  production  can  be  estimated  through  the 
Obukhov  length  (see  Results). 


Table  1.  Summary  of  tracer  release  data. 


Trial 

Date 

Start  Time 
(MDT) 

Q 

(1  min“^) 

(m,m,m) 

1 

21-09-2001 

01:04:00 

175 

(24,19,10,05,0.15) 

2 

21-09-2001 

02:51:00 

200 

(24.19,-30.15,0.15) 

3 

24-09-2001 

18:52:00 

200 

(39.51,90.45,0.15) 

4 

24-09-2001 

19:35:00 

200 

(39.51,90.45,1.8) 

5 

24-09-2001 

20:34:00 

200 

(24.19,50.25,1.8) 

6 

24-09-2001 

21:02:00 

200 

(24.19,30.15,0.15) 

7 

24-09-2001 

21:51:00 

200 

(29.64,-30.15,0,15) 

8 

24-09-2001 

22:13:00 

200 

(29.64,-30.15,1.8) 

9 

24-09-2001 

22:35:00 

200 

(31.85,-30.15,2.6) 

10 

24-09-2001 

23:04:00 

200 

(8.87,0.0,1.8) 

11 

25-09-2001 

18:30:00 

225 

(8.87, 70.35,1.8) 

12 

25-09-2001 

18:49:00 

225 

(8.87, 70.35,0.15) 

13 

25-09-2001 

21:51:00 

225 

(14.32,70.35,1.8) 

14 

25-09-2001 

22:56:00 

225 

(29.64,90.45,0.15) 

15 

25-09-2001 

23:21:00 

225 

(31.85,90.45,2.6) 

16 

25-09-2001 

23:54:00 

225 

(31.85,50.25,5.2) 

17 

26-09-2001 

20:55:00 

225 

(-24.0,-70.35,1.3) 

18 

26-09-2001 

21:32:00 

225 

(-1.0,-90.45,1.3) 

19 

26-09-2001 

21:58:00 

225 

(1,21.-90.34,2.6) 

20 

26-09-2001 

22:24:00 

225 

(-1.0,-70.35,1.3) 

21 

26-09-2001 

22:51:00 

225 

(-24,0,-70.35,1.3) 

The  pre-processing  of  the  concentration  time  series  consisted  of  two  basic  steps.  (1)  drift 
correction  to  the  baseline  level  of  the  concentration  time  series  to  remove  any  slow  variations 
(trends)  in  the  level  over  the  sampling  period;  and,  (2)  conversion  of  the  concentration  time 
series  (in  A-D  units)  to  absolute  concentration  (in  ppm  by  volume)  by  application  of  the 
appropriate  calibration  curve  (modelled  as  a  second-degree  polynomial)  to  the  data.  The 
baseline  drift  was  removed  by  least-squares  fitting  cubic  splines  to  the  mean  levels  of  the 
noise  baseline  segments  (segments  containing  no  rapid  fluctuations  of  significantly  non-zero 
concentration)  in  the  concentration  time  series.  The  concentration  detectors  were  calibrated 


10 


DRDC  Sufficld  TR  2003-097 


at  regular  intervals  (throughout  the  experimental  period)  over  their  entire  operating  range 
using  reference  concentrations  of  0.  5.  10.2,  50,  110,  200,  and  400  ppm  propylene  in  air.  The 
second-order  polynomial  used  to  characterize  the  calibration  curve  was  found  to  provide  an 
accurate  representation  (to  within  ±2  %)  for  the  calibration  data  over  the  entire  operating 
range  of  the  detectors. 

Almost  all  the  experiments  in  this  study  were  conducted  in  a  slightly  to  very  stably  strat¬ 
ified  atmospheric  surface  layer.  For  stable  conditions,  the  crosswind  velocity  variance  can 
increase  substantially  due  to  the  decrease  of  the  friction  velocity  and  increase  of  the  relative 
contributions  from  the  mesoscale  variations.  This  can  result  in  larger  scale  two-dimensional 
motions  (e.g..  larger  than  the  characteristic  length  or  width  of  the  obstacle  array)  that  may 
persist  for  long  periods  in  the  stable  boundary  layer,  and  these  motions  can  enhance  the 
meandering  motion  of  plumes  (Hanna  [27];  Etling  [28]).  Hence,  for  long  sampling  times, 
these  larger  two-dimensional  motions  will  dominate  the  dispersion  obscuring  the  effect  of 
the  obstacle  array  on  the  turbulent  diffusion.  Because  the  crosswind  velocity  variance  is 
sensitive  to  the  range  of  scales  included  in  its  determination,  the  measurement  of  plume 
spread  under  these  conditions  can  be  potentially  very  sensitive  to  the  choice  of  sampling 
time. 

For  a  sampling  time  of  approximately  15  min  used  in  our  experiments,  the  plume  meander¬ 
ing  caused  by  the  quasi-horizontal  mesoscale  eddies  in  the  stable  boundary  layer  resulted 
generally  in  a  slow  drift  in  direction  of  the  plume  across  the  horizontal  sampling  lines.  This 
result  is  generally  further  exacerbated  by  the  difficulties  and  complexities  associated  with 
unsteadiness,  non-uniformity,  and  sensitivity  of  flows  to  terrain  slope  in  the  stable  boundary 
layer.  To  remove  the  effects  of  this  non-stationarity  from  the  data,  the  concentration  time 
series  for  each  plume  dispersion  experiment  were  conditionally  sampled  based  on  a  careful 
examination  of  the  records  of  mean  wind  speed  and  direction  provided  by  the  2-D  sonic 
anemometers  deployed  on  the  upwind  mast.  These  records  of  mean  wind  speed  and  direc¬ 
tion  were  used  to  extract  the  200  s  period  within  each  15  min  plume  dispersion  experiment 
that  exhibited  the  least  variation  in  the  mean  wind  speed  and  direction.  This  conditional 
sampling  of  the  velocity  and  concentration  time  series  restricts  the  range  of  scales  included 
in  its  determination,  and  also  makes  possible  a  meaningful  comparison  of  our  results  with 
previous  investigations,  both  for  field  and  wind  tunnel  experiments  (e.g.,  Macdonald  et  al. 
[12]). 

Our  experimental  measurements  correspond  to  near-field  dispersion  in  the  obstacle  array, 
where  the  travel  time  T  is  generally  much  less  than  the  finite  duration  r  (?»  15  min)  of 
the  release  (i.e..  T  ^  r).  In  all  the  selected  experiments,  T  <  Tg  <C  r  where  Tg  =  200  s 
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is  the  fixed  sampling  duration  of  the  standard  period  exhibiting  the  least  variation  in  the 
mean  wind.  For  the  short  range  of  downwind  distances  covered  in  our  experiments,  we 
expect  the  effective  \'ariation  in  the  turbulent  velocity  statistics  that  determine  the  plume 
dispersion  to  be  dominated  by  the  fluctuations  at  the  point  of  release  over  the  period  of  the 
duration  T,  of  sampling  of  the  concentration  field,  owing  to  the  fact  that  the  Lagrangian 
correlation  coefficient  is  expected  to  fall  only  slowly  (Hay  and  Pasquill  [29]).  In  other 
words,  the  magnitude  of  the  scale  of  turbulence  is  expected  to  be  irrelevant  to  the  turbulent 
spread  of  contaminant  material  close  to  a  continuous  point  source.  However,  the  scale 
of  the  turbulence  must  become  increasingly  important  to  the  determination  of  the  plume 
spreading  as  the  ratio  of  the  time  of  travel  to  the  Lagrangian  time  scale  increases. 

It  seems  reasonable  that  the  turbulence  statistics  (e.g.,  velocity  variances)  that  determine 
the  plume  dispersion  over  the  range  of  limited  downwind  distances  covered  in  our  exper¬ 
iments  will  be  those  corresponding  to  the  sampling  period  Tg  of  the  plume  concentration 
field,  provided  the  time  of  travel  T  is  not  substantially  larger  than  Tg.  In  addition,  in 
some  of  the  experiments  (especially  those  conducted  under  near  neutral  or  slightly  stable 
stratification  with  the  higher  mean  wind  speeds),  we  found  a  remarkable  steadiness  in  the 
statistical  values  of  the  standard  deviation  of  velocity  (e.g.,  in  these  experiments,  we  found 
that  the  crosswind  velocity  standard  deviation  sampled  over  600  s  was  less  than  15  %  larger 
than  that  sampled  over  200  s,  implying  only  a  very  gradual  increase  in  the  velocity  standard 
deviation  with  increasing  sampling  time). 
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Results 


Velocity  statistics  upwind,  within,  and  downwind  of  the  array 


A  compendium  of  mean  wind  and  turbulence  statistics  data  for  the  sampling  periods  selected 
for  each  experiment  is  given  in  Table  2.  Here,  wind  components  are  calculated  in  the  rotated 
coordinate  system  with  streamwise  (u),  spanwise  (u),  and  vertical  {w)  components  aligned 
with  the  X-,  y-  and  ^-directions,  respectively  (and  their  standard  deviations  are  given  by 
(Tu,  ay,  and  a^,  respectively).  Measurements  of  the  mean  horizontal  wind  speed  Sq4,  mean 
wind  direction  a,  and  standard  deviation  in  the  mean  wind  direction  a^  obtained  from  the 

2- D  sonic  anemometer  mounted  at  a  height  of  4  m  on  the  upwind  mast  are  summarized 

in  Table  2.  Here,  S  =  (u'^  -f  is  the  horizontal  wind  speed  where  the  overbar  is  used 

to  denote  a  time  averaged  value,  and  a  =  tan“^(u/u)  is  the  mean  wind  direction  with  a 
positive  angle  measured  anti-clockwise  from  the  horizontal  o’-axis  (which  is  aligned  along 
the  direction  normal  to  the  front  face  of  the  array).  Hence,  a  =  0°  corresponds  to  a  mean 
wind  direction  normal  to  the  obstacle  array.  Measurements  of  turbulence  obtained  with  the 

3- D  sonic  anemometer/thermometer  at  the  4-m  level  on  the  32-m  vertical  tower  near  the 
array  centre  are  summarized  also  in  Table  2.  The  turbulence  kinetic  energy  (TKE),  k,  was 
estimated  from 

k  =  +  (tI  +  al) .  (1) 

The  Obukhov  length,  Lmo-  was  estimated  as 


Amo  = 


ulT 

ngw'T'  ’ 


(2) 


where  u*  is  the  friction  velocity,  T  is  the  mean  sonic  temperature,  k  ^  0.4  is  the  von 
Karman  constant  for  momentum  diffusion,  g  is  the  acceleration  due  to  gravity,  and  w'T'  is 
the  vertical  temperature  flux.  The  influence  of  the  mean  humidity  flux  on  the  value  of  Amo 
was  considered  to  be  negligible.  The  friction  velocity  was  estimated  from 


u* 


+  {y'w'Y 


-,1/4 


(3) 


Finally,  the  mean  rate  of  dissipation  of  turbulence  kinetic  energy,  e.  was  estimated  from  the 
measured  u*  and  Amo  using  similarity  theory  modified  for  use  in  stable  conditions  as 


/I  ,  4  \ 

K  \ z  Amo / 


(4) 
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Note  that  Eq.  (4)  was  obtained  from  the  TKE  transport  equation  assuming  local  equilibrium 
whereby  a  balance  exists  between  (shear  and  buoyant)  production  and  dissipation  of  TKE 
{uIdU /dz  -  ui/inL-isio)  =  c.  where  U  is  mean  wind  speed),  with  the  mean  wind  shear  for 
stable  conditions  given  by 

where  (3  (Dyer  [30]). 


Table  2.  Summary  of  velocAty  statistics. 


Trial 

So4 

a 

(Jet 

k 

Lmo 

e 

(ms~') 

(dag) 

(deg) 

(m^s~^) 

(ms“') 

(m) 

(m^s“^) 

1 

2.35 

17 

5.7 

0.359 

0.26 

91 

0.013 

2 

2.01 

30 

7.3 

0.306 

0.25 

62 

0.013 

3 

3.06 

-49 

8.6 

0.436 

0.32 

330 

0.022 

4 

1.63 

-48 

4.0 

0.148 

0.08 

5.8 

0.0013 

5 

2.69 

-26 

1.8 

0.251 

0.17 

4.8 

0.012 

6 

1.89 

-10 

3.8 

0.218 

0.16 

7.7 

0.0073 

7 

2.30 

36 

9.8 

0.409 

0.35 

150 

0.030 

8 

2.68 

30 

7.9 

0.428 

0.35 

150 

0.030 

9 

10 

2.32 

2.56 

36 

17 

8.2 

6.2 

0.387 

0.367 

0.26 

0.25 

48 

74 

0.014 

0.012 

11 

7.93 

-41 

9.5 

1.46 

1.1 

28000 

0.80 

12 

7.26 

-50 

8.5 

0.877 

0.76 

2500 

0.28 

13 

3.89 

-41 

7.1 

0.402 

0.46 

140 

0.066 

14 

5.02 

-42 

9.2 

0.877 

0.66 

240 

0.19 

15 

4.55 

-39 

8.5 

0.718 

0.50 

170 

0.087 

16 

4.49 

-47 

7.9 

0.727 

0.44 

120 

0.059 

17 

3.34 

39 

8.3 

0.362 

0.36 

170 

0.032 

18 

4.00 

39 

7.6 

0.582 

0.42 

220 

0.048 

19 

2.98 

43 

9.1 

0.505 

0.39 

130 

0.042 

20 

2.63 

26 

9.5 

0.484 

0.35 

120 

0.030 

21 

3.38 

36 

7.8 

0.537 

0.37 

130 

0.036 

Mean  velocities  and  turbulence  statistics  were  recorded  by  various  sonic  anemometers  up¬ 
wind,  through,  and  downwind  o.^  the  obstacle  array.  The  sonic  anemometer  on  the  upwind 
mast  at  a.  height  of  4  m  wa.^^  used  to  provide  the  reference  mean  horizontal  wind  speed 
504  and  turbulence  kinetic  energy  ko4  for  normalizing  the  mean  velocity  and  turbulence 
statistics,  respectively,  measured  within  and  downwind  of  the  obstacle  array.  Note  that  the 
subscript  “04”  is  used  here  to  denote  a  flow  cjuantity  that  is  measured  in  the  undisturbed 
upwind  flow  at  the  4-ni  level.  Because  the  upwind  2-axis  sonic  anemometer  at  the  4-m 
level  provided  only  the  two  horizontal  components  of  the  velocity,  the  reference  turbulence 
kinetic  energy  koA  was  estimated  from  the  measurements  of  and  a„,  and  an  approxima¬ 
tion  that  (T,,.  ~  0.4(cr^  +  The  estimate  for  cTu,.  was  obtained  from  using  the  usual 
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turbulence  scaling  ratios  typical  of  a  neutrally  stratified  surface  layer;  namely,  =  2.5, 

Gy/Ut:  =  2.0,  and  =  1.25  (u,  is  the  friction  velocity). 

The  spatial  variation  of  the  normalized  mean  horizontal  wind  speed  S/Sq^  at  various  heights 
upwind,  within,  and  downwind  of  the  obstacle  array  is  exhibited  in  Figure  2  for  a  mean  wind 
direction  nearly  perpendicular  (|q|  =  10°)  and  at  a  large  flow  incidence  angle  (|a|  =  41°) 
to  the  front  face  of  the  array.  The  reduction  in  the  mean  wind  speed  below  or  at  canopy 
height  is  more  significant  for  the  near  normally-incident  flow  than  for  the  flow  at  an  oblique 
incidence,  although  in  both  cases  minimum  wind  speeds  at  or  below  the  canopy  top  were 
observed  after  the  third  row.  However,  above  the  canopy  height,  the  situation  is  reversed 
with  the  reduction  in  the  mean  wind  speed  being  more  pronounced  for  the  41°  wind  direction 
than  for  the  near  perpendicular  wind  direction.  In  addition,  the  mean  wind  shear  observed 
near  the  canopy  top  in  the  array  for  a  nearly  perpendicular  wind  was  larger  than  that 
observed  for  an  oblique  flow  incidence. 

Figure  3  shows  the  spatial  variation  of  the  normalized  turbulence  kinetic  energy  k/ko^  with 
depth  into  the  obstacle  array  at  various  heights  for  winds  at  near  perpendicular  and  at  an 
oblique  flow  incidence.  The  TKE  increases  at  or  near  the  canopy  top  as  one  progresses 
through  the  obstacle  array.  This  is  due  to  an  intense  shear  layer  at  the  top  of  the  canopy, 
where  the  kinetic  energy  of  the  mean  flow  is  converted  into  turbulence  kinetic  energy  (shear 
production  of  TKE).  In  addition,  the  TKE  generated  upwind  of  a  given  fixed  point  in 
the  array  will  be  exported  downwind  by  local  advection  and  turbulent  (and/or  pressure) 
transport.  We  note  that  the  increase  in  the  TKE  through  the  obstacle  array  (below,  at, 
and  just  above  the  canopy  top)  is  considerably  greater  for  the  case  of  near  normal  flow 
incidence  than  for  that  of  an  oblique  flow  incidence.  The  variation  in  the  TKE  below 
and  above  the  canopy  height  is  significantly  greater  for  a  near  normal  wind  than  for  an 
oblique  wind  direction.  Consequently,  it  would  appear  that  the  turbulence  generated  either 
by  the  obstacles  themselves  in  the  wake  flow  (wake  production  of  TKE)  or  by  the  shear¬ 
generated  turbulence  produced  near  the  canopy  top  in  the  inflected  mean  velocity  profile 
is  considerably  greater  for  winds  perpendicular  to  rather  than  at  some  oblique  angle  to  the 
front  of  the  obstacle  array. 

Figure  4  displays  the  variation  of  the  bulk  drag  coefficient  (=  u^/SD  as  a  function 
of  stability  parameter  using  sonic  anemometer  data  acquired  at  a  height  of  4  m 

on  the  32-m  vertical  tower  located  near  the  array  centre.  For  L^iq  ^  0.01,  the  drag 
coefficient  diminishes  with  increasing  stability  (decreasing  L^io)-  Note  that  the  stability 
regime  >  0.01  corresponds  to  a  very  stable  stratification  where  there  is  significant 

suppression  of  mechanical  turbulence  by  the  stable  stratification  (at  4-m  height  near  the 
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array  centre).  As  increases  within  this  range  of  values,  the  turbulence  in  the  stable 

surface  layer  (overlying  the  obstacle  array)  becomes  more  intermittent  and  patchy,  allowing 
the  upper  portions  of  the  layer  to  decouple  from  the  surface  forcings.  In  particular,  for 
-^ivio  ~  0-01  coefficient  exhibits  a  power-law  decay  with  increasing  stability  (Cd  ~ 

(L^jo)"’')  where  the  power-law  exponent  is  found  to  be  r  «  1/2  (±10  %)■ 

The  mean  wind  speed  profile,  adjusted  using  the  log-linear  profile  for  stable  conditions,  can 
be  written 

where  Z  =  z  —  d  is  the  displaced  (effective)  height,  zq  is  the  aerodynamic  roughness  length, 
d  is  the  displacement  height,  and  u,  is  the  surface  friction  velocity.  The  parameter  (3  is 
usually  taken  to  be  5,  although  published  values  can  range  from  3  to  9  (Businger  et  al.  [31]; 
Dyer  [30];  Hogstrdm  [32]).  Figure  4  shows  the  drag  coefficient  predicted  using  Equation  (5) 
with  /3  =  5  and  the  upwind  roughness  length  and  displacement  height  zq  —  0.045  m  and 
d  =  0.37  m  (see  Site  and  Topography).  Note  that  the  value  of  the  drag  coefficient  is 
underestimated  using  the  upwind  roughness  length.  This  is  because  the  air  flow  above  the 
obstacle  array  produces  an  additional  momentum  flux,  the  form  (pressure)  drag,  resulting 
in  an  increased  friction  velocity,  as  =  (T^±r^)V2,  ^^ere  r  is  the  total  kinematic  shear 
stress,  is  the  turbulent  kinematic  shear  stress,  and  is  the  additional  kinematic  shear 
stress  arising  from  the  form  drag  due  to  the  presence  of  the  obstacles.  At  a  given  height  and 
wind  speed  and  assuming  that  the  mean  wind  profile  above  the  obstacle  array  retains  the 
form  given  in  Equation  (5).  this  implies  an  increased  roughness  parameter  for  the  obstacle 
array  {zq  increases  with  increasing  Cd)-  In  particular,  the  bulk  drag  coefficient  appears 
to  be  consistent  with  Equation  (5)  with  3  =  5  and  an  increased  aerodynamic  roughness 
length  zq  =  0.20  m.  For  the  obstacle  arra\-.  the  displacement  height  is  assumed  to  be  70% 
of  the  obstacle  height  H,  so  d  ^  1.8  m  (reflecting  an  increased  displacement  of  the  flow 
over  the  obstacle  array).  The  increased  aerodynamic  roughness  length  here  simply  reflects 
the  two  contributions  to  the  total  stress;  namely,  the  turbulent  surface  friction  (which  can 
be  described  by  the  roughness  length  of  the  upwind  fetch)  and  the  pressure  or  form  drag 
(which  is  produced  by  the  obstacle  array). 


Mean  concentration 


Figure  5  displays  typical  crosswind  profiles  of  the  mean  concentration  C  in  the  obstacle 
array  at  two  downwind  distances  xl  from  the  source  for  near  normal  flow  incidence  (TYial 
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6).  Here,  xl  refers  to  the  Euclidean  distance  in  the  horizontal  x-y  plane  between  the  point 
source  location  and  the  position  of  the  plume  centreline  along  the  horizontal  sampling 
line.  Also,  in  Figure  5,  a  spanwise  coordinate  Y  referred  to  a  survey  line  parallel  to  the 
lengthwise  (or.  along  the  x-direction)  edge  of  the  array  and  4.76  m  west  of  the  obstacles 
in  column  8  has  been  used;  thus,  Y  =  y  +  81.21  m.  The  measurements  suggest  that  the 
crosswind  concentration  profiles  are  well  represented  by  a  Gaussian  distribution,  a  result 
that  is  consistent  with  the  findings  of  Davidson  et  al.  [9]  and  Macdonald  et  al.  [11]. 

Vertical  profiles  of  the  mean  concentration  at  the  plume  centreline  are  shown  in  Figure  6  for 
two  different  downwind  distances  and  source  release  heights  for  array  plumes  from  Trials  6 
and  10.  Depending  on  the  location  and  height  of  the  source  and  on  the  downwind  distance 
from  the  source,  the  maximum  mean  concentration  in  a  vertical  plume  cross-section  can 
occur  at  the  surface  [Figure  6(a)]  or  above  the  surface  [Figure  6(b)].  Hence,  vertical  profiles 
of  mean  concentration  for  plumes  dispersing  through  an  obstacle  array  develop  in  a  rather 
complex  manner  in  contrast  to  the  simple  Gaussian  crosswind  profiles. 

To  characterize  the  lateral  and  vertical  growth  of  the  mean  plume,  we  fitted  a  Gaussian 
profile  of  the  form 

C(!/)=Co,exp(^-fc^),  (6) 

to  the  lateral  profiles  of  mean  concentration  and  a  reflected  Gaussian  model  of  the  form 

CW  =  C„,|exp(-<i^)+cxp(-ii^)},  (7) 

to  the  vertical  profiles  of  mean  concentration.  Here,  Coy  and  Coz  are  the  lateral  plume 
centreline  mean  concentration  and  vertical  plume  centreline  mean  concentration  (neglect¬ 
ing  reflection  from  the  ground),  respectively;  Uy  and  Uz  are  the  lateral  and  vertical  plume 
standard  deviations,  respectively:  and,  yc  and  Zc  are  the  lateral  displacement  of  the  plume 
centroid  and  the  plume  height  centroid,  respectively.  Figures  5  and  6  show  typical  examples 
of  best-fit  Gaussian  and  reflected  Gaussian  forms  to  the  lateral  and  vertical  mean  concen¬ 
tration  profiles,  respectively.  Although  the  Gaussian  form  provides  very  good  agreement 
with  the  measured  lateral  mean  concentration  profiles,  the  reflected  Gaussian  form  provides 
a  less  than  ideal  description  of  the  vertical  mean  concentration  profiles  for  many  cases.  In 
consequence,  the  derived  vertical  plume  spread  using  the  reflected  Gaussian  model,  while 
useful  for  both  modelling  and  practical  purposes,  is  not  as  accurate  as  the  derived  lateral 
plume  spread  using  the  Gaussian  model. 

Releases  within  or  just  upwind  (in  the  displacement  zone  just  upwind  of  the  front  row)  of 
the  obstacle  array  can  result  in  significantly  different  vertical  mean  concentration  profiles 
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through  the  plume  than  releases  well  upwind  of  the  array  (well  upwind  of  the  zone  where  the 
flow  is  first  influenced  by  the  presence  of  the  array).  Figure  7  shows  vertical  profiles  of  the 
non-dimensional  mean  concentration  K  (=  10“®C'5o4-f^^/Q;  vhere  C  is  mean  concentration 
in  ppm  (i.e.,  for  pure  undiluted  gas  10~^C  =  1),  5o4  is  the  upwind  horizontal  mean  wind 
speed  at  4  m,  and  Q  is  the  tracer  volume  release  rate  in  m^  for  these  two  release 
configurations.  The  Figure  clearly  illustrates  the  differences  in  the  mean  plume  dispersion 
for  these  two  situations  (taken  from  Trials  17  and  18).  Both  releases  were  for  a  source  set 
at  a  height  of  H/2  with  the  incident  wind  direction  at  a  =  39°.  For  releases  well  upwind  of 
the  front  of  the  array,  the  mean  flow  streamlines  are  accelerated  and  deflected  over  the  top 
of  the  array  causing  the  plume  centroid  to  be  elevated  and  increasing  the  vertical  growth  of 
the  plume.  In  contrast,  for  releases  within  or  just  upwind  of  the  array,  the  downwash  of  the 
plume  material  around  the  obstacles  arising  from  the  diffusion  of  vorticity  from  the  wakes 
of  the  obstacles  results  in  an  observed  maximum  mean  concentration  at  ground  level.  The 
effective  lifting  of  the  plume  centroid  for  sources  located  upwind  of  an  obstacle  array  has 
previously  been  reported  by  Davidson  et  al.  [9]. 

The  presence  of  the  obstacles  results  in  a  lateral  deflection  in  the  direction  of  the  mean 
plume  centreline  relative  to  the  mean  wind  direction  in  the  approach  flow.  Figure  8,  which 
shows  the  mean  plume  centreline  angle  6  (the  angle  between  the  positive  x-direction  and  the 
vector  from  the  source  location  to  the  location  of  the  mean  plume  centreline)  as  a  function 
of  the  mean  wind  direction  a  in  the  approach  flow,  illustrates  this  phenomenon.  Note  that 
for  near-normal  flow  incidence  with  |q1  <  20°,  the  plume  centreline  direction  is  deflected 
towards  the  normal  to  the  front  face  of  the  array  (|0|  <  |o;|).  This  appears  to  be  due  to  the 
wind  deflection  or  channeling  down  the  street  canyons  of  the  obstacle  array  aligned  roughly 
along  the  SE-NW  direction  (i.e.,  more  precisely  the  wind  direction  within  the  array  should 
be  around  30°  from  true  north).  However,  for  a  greater  oblicjuity  of  flow  incidence,  the  mean 
plume  centreline  direction  is  deflected  away  from  the  normal  to  the  front  face  of  the  array, 
implying  the  plume  lateral  displacement  is  due  to  the  channeling  of  the  wind  along  the 
street  canyons  of  the  array  that  are  aligned  approximately  along  the  SW-NE  direction  (i.e., 
more  precisely  the  wind  direction  within  the  array  should  be  around  60  or  120°  from  true 
north).  In  summary  then,  the  surface  flow  tends  to  be  channeled  approximately  parallel  to 
the  street  canyons  roughly  in  the  SE-N^^  direction  when  the  wind  above  the  obstacle  array 
has  a  large  component  in  this  direction  (when  the  approach  flow  is  approximately  normal  to 
the  front  face  of  the  array),  but  when  the  wind  aloft  of  the  array  has  a  significant  component 
away  from  this  direction  (when  the  approach  flow  is  at  an  oblique  angle  greater  than  about 
20°  to  the  array  front)  then  the  wind  in  the  array  tends  to  be  channeled  approximately 
parallel  to  the  street  canyons  roughly  in  the  SW-NE  direction  (parallel  to  the  front  face 
of  the  array).  This  channeling  of  the  surface  wind  in  the  obstacle  array  is  reflected  in  the 
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observed  behaviour  of  the  lateral  displacements  in  the  mean  centreline  directions  of  the 
array  plume. 

The  measured  rates  of  lateral  and  vertical  plume  growth  in  the  obstacle  array  for  near 
perpendicular  flow  incidence  are  illustrated  in  Figure  9.  Here,  it  should  be  noted  that 
the  vertical  plume  growth  is  represented  by  the  vertical  extent,  az^t  =  (^z  +  ^c)  ^  )• 
To  provide  a  reference  for  comparison,  conventional  results  for  lateral  and  vertical  plume 
growth  in  open  terrain  obtained  for  various  Pasquill-Gifford  stability  classes  B,  C,  D,  and 
E  following  Briggs  (for  a  10-minute  sampling  duration)  are  also  included  (Stern  et  al.  [33]). 
Note  that  the  lateral  and  vertical  plume  spreads  as  well  as  downwind  distance  from  the 
source  for  the  unobstructed  plume  (where  the  only  scale  factors  are  roughness  length  and 
boundary-layer  height)  have  been  normalized  by  FT  in  order  to  allow  a  comparison  to  be 
made  with  the  array  plume  results.  At  the  short  ranges  considered  here  where  the  statistical 
theory  (Csanady  [34])  predicts  that  lateral  and  vertical  plume  spreads  are  proportional  to 
downwind  distance  from  the  source,  it  is  expected  that  the  scaling  of  these  plume  parameters 
should  be  similar  over  a  broad  range  of  scales  (providing  a  justification  for  the  use  of  if  as 
the  scaling  parameter  for  the  unobstructed  plume  parameters). 

The  measurements  of  plume  dispersion  in  the  array  shown  in  Figure  9  were  conducted  under 
moderately  stable  background  atmospheric  conditions  (stability  class  E),  so  the  lateral  and 
vertical  plume  spreads  for  the  array  plume  should  be  compared  with  the  open-terrain  Briggs 
dispersion  curve  for  stability  class  E.  This  comparison  shows  that  the  lateral  and  vertical 
plume  spreads  of  the  array  plume  are  both  increased  relative  to  that  of  the  unobstructed 
plume  dispersion  under  the  same  atmospheric  stratification.  At  x^/H  ks  10,  the  lateral 
plume  spread  is  similar  in  magnitude  to  that  in  a  plume  dispersing  in  open  terrain  under 
stability  class  B  and  the  subsequent  rate  of  lateral  plume  growth  in  the  array  is  similar 
to  that  of  an  open-terrain  plume  in  stability  class  D.  Similarly,  the  vertical  plume  spread 
through  the  obstacle  array  at  Xl/H  «  15  is  similar  to  that  of  an  open-terrain  plume  in 
stability  class  C,  although  the  subsequent  rate  of  growth  is  comparable  to  that  for  stability 
class  E  (viz.,  approximately  the  same  as  that  for  an  unobstructed  plume  dispersing  under  the 
same  background  atmospheric  stratification).  Qualitatively  similar  results  were  obtained 
by  Macdonald  et  al.  [12]  for  increases  in  both  the  lateral  and  vertical  plume  widths  for  wide 
obstacles  with  aspect  ratios  W/H  greater  than  unity.  There,  as  here,  it  appears  that  for 
wide  obstacles  with  W/H  >  1  up  to  two-dimensional  obstacles  {W/H  — >  oo),  the  lateral 
plume  spread  increased  by  a  factor  of  2-4  (at  leeist  near  the  source  in  the  “inflationary 
phase”  of  plume  growth). 

A  number  of  physical  mechanisms  are  responsible  for  the  observed  rapid  initial  increase 
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(see  Figure  9)  in  the  lateral  and  vertical  spreads  of  the  array  plume  ("inflationary  phase  ’  of 
plume  growth).  The  lateral  divergence  of  streamlines  around  wide  and  low  obstacles  followed 
by  entrainment  and  rapid  mixing  of  the  plume  material  across  the  width  W  of  the  obstacle 
in  the  wake  region  and/or  the  lateral  plume  displacement  (and  the  concomitant  increased 
lateral  shear  effect)  increases  the  lateral  plume  spread.  On  the  other  hand,  the  rapid  mixing 
of  plume  material  in  the  obstacle  wakes  up  to  the  height  H  of  the  obstacles  and  the  vertical 
divergence  of  mean  flow  streamlines  near  the  front  of  the  array  increases  the  vertical  plume 
spread.  In  the  case  of  a  release  upwind  of  the  front  face  of  the  array,  the  upward  drift 
(deflection)  of  the  mean  flow  streamlines  near  the  leading  edge  of  the  array  results  in  an 
effective  lifting  of  the  plume  centroid  which  enhances  the  upward  vertical  diffusion.  In  turn, 
plume  material  swept  upwards  to  the  obstacle  height  by  the  wake  motions  can  sometimes  be 
subsequently  entrained  into  intermittent,  large  eddy  structures  that  can  rapidly  transport 
the  material  upwards  to  a  signiflcant  height.  It  is  conceivable  that  these  canopy-scale 
eddies  dominate  the  vertical  transport  of  momentum  and  scalar  in  a  canopy.  These  large- 
scale  eddy  structures  are  the  result  of  the  inviscid  instability  to  small  perturbations  arising 
from  the  inflected  canopy  mean  velocity  profile  which  sets  the  pattern  for  the  coherent  eddy 
generation  at  the  canopy  top.  The  characteristics  of  these  coherent  eddies  have  been  studied 
by  Raupach  et  al.  [35]  using  the  canopy-mixing  layer  analogy,  and  have  been  shown  to  scale 
with  H  rather  than  z. 


The  observed  behaviour  of  the  array  plume  growth  rate  can  be  explained  as  follows.  The 
rate  of  spread  of  the  plume  in  the  lateral  and  vertical  directions  may  be  expressed  in  terms 
of  the  lateral  (A'y)  and  vertical  [Kz)  turbulent  diffusivities,  respectively,  as 


K,i  j  doz  _  Kz 

dx  Say  dx  Suz 


(8) 


where  Ky  «  u.™sA  and  Kz  ^  (u,.,,s  and  uv„,s  are  characteristic  velocities  of  the 

turbulence  in  the  lateral  and  vertical  directions,  respectively;  and.  A  is  the  turbulence 
integral  length  scale).  The  turbulence  intensity  (urms/S'  or  w^ms/S)  is  expected  to  be  larger 
in  the  "urban"  canopy  than  in  open  terrain  due  to  the  obstacle-generated  turbulence.  In 
order  for  the  array  plume  growth  rate  to  be  comparable  to  that  in  an  open-terrain  plume, 
the  increase  of  the  turbulence  intensity  in  the  “urban”  canopy  relative  to  that  in  the  open 
terrain  (which  increases  the  plume  growth  rate)  must  be  compensated  by  the  effects  of  the 
following  two  factors.  Firstly,  the  scale  of  the  turbulence  in  the  canopy  must  be  reduced 
relative  to  that  in  the  open  terrain.  Secondly,  the  array  plume  spread  (cr.y  or  cr^)  must 
be  increased  relative  to  the  open-terrain  plume  spread.  Both  of  these  factors  decrease  the 
plume  growth  rate.  In  particular,  the  turbulence  length  scale  in  the  canopy  is  comparable 
to  the  height  of  the  obstacles,  whereas  the  turbulence  length  scale  in  the  open  terrain  (m 
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the  incident  flow)  is  generally  much  larger  than  typical  building  (obstacle)  dimensions.  For 
the  obstacle  array  studied  here,  these  three  effects  (increased  turbulence  intensity  on  the 
one  hand,  and  decreased  turbulence  integral  length  scale  and  increased  plume  spreads  on 
the  other  hand)  compensate  each  other  with  the  result  that  the  rate  of  plume  growth  in 
the  array  after  the  initial  rapid  increase  (“inflationary  phase")  is  comparable  to  that  in  an 
open-terrain  plume  (at  least  over  the  range  of  xl/H  for  which  we  have  measurements). 
However,  given  that  doy/dx  oc  (7“^  [Equation  (8)],  it  is  conceivable  that  at  very  large 
downstream  distance  {xl/H  »  1)  the  lateral  plume  spread  rate  of  the  array  plume  could 
be  actually  smaller  than  that  of  the  unobstructed  plume  (owing  to  the  larger  initial  lateral 
array  plume  spread).  This  would  imply  that  the  lateral  spread  of  the  open-terrain  plume 
would  eventually  “catch  up”  to  that  of  the  array  plume,  so  that  at  this  stage  of  plume 
development  the  open-terrain  and  array  plumes  would  be  comparable  in  “size” . 

Figure  10  exhibits  the  decay  of  the  non-dimensional  mean  concentration  at  the  plume  cen¬ 
treline  Kq  with  downwind  distance  from  the  source  through  the  obstacle  array  for  different 
source  heights  and  various  directions  6  of  plume  centreline  dispersion  through  the  array. 
The  open  terrain  measurements  of  Kq  shown  here  were  obtained  by  combining  the  results 
from  a  previous  series  of  field  experiments  (Yee  et  al.  [24];  Yee  et  al.  [25])  conducted  under 
near-neutral  stability  conditions.  Note  that  this  measured  variation  of  Kq  in  open  terrain 
under  near-neutral  atmospheric  stratification  is  generally  consistent  with  the  variation  of 
the  Pasquill-Gifford  neutral  (D)  curve.  However,  the  non-dimensional  plume  centreline 
mean  concentration  in  the  obstacle  array  (obtained  under  moderately  to  very  stable  strati¬ 
fication)  falls  nearer  the  stability  class  C  curve  rather  than  the  class  E  curve.  Hence,  plume 
centreline  mean  concentrations  are  reduced  considerably  in  the  obstacle  array  relative  to 
the  open-terrain  plume  results.  This  effect  for  the  wide  obstacle  array  [WfH  «  4.8)  studied 
here  is  different  than  that  obtained  for  an  array  of  cubes  {W/H  —  1)  where  it  was  found 
that  the  plume  centreline  mean  concentration  in  the  array  is  similar  to  the  open-terrain 
plumes  (Davidson  et  al.  [9];  Macdonald  et  al.  [11]).  It  is  evident  that  for  the  wide  ob¬ 
stacle  array  considered  here,  the  reduction  in  the  mean  wind  speed  within  the  array  does 
not  compensate  for  the  significant  increases  in  both  the  lateral  and  vertical  plume  widths 
(and  significant  increase  in  the  plume  area)  and  therefore  the  centreline  mean  concentration 
Co  in  the  array  plume  is  significantly  smaller  than  the  open-terrain  plume  under  similar 
atmospheric  conditions. 

It  is  interesting  to  note  that  at  downwind  distances  from  the  source  greater  than  about 
Xl/H  ^  30,  the  rate  of  decrease  of  the  centreline  mean  concentration  in  the  array  plume 
is  less  than  that  of  an  open-terrain  plume  dispersing  under  similar  atmospheric  conditions. 
This  can  be  understood  as  follows.  For  simplicity,  assume  that  the  mean  concentration  field 
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can  be  approximated  by  a  Gaussian  distribution  in  both  the  lateral  and  vertical  directions 
(i.e.,  neglect  the  surface  reflection),  with  a  centerline  mean  concentration  given  by  Cq  = 
Q^/{2-KSaya^),  where  Qm  is  the  quantity  of  matter  released  per  unit  time  (Stern  et  al. 
[33]).  Consequently,  the  rate  of  decrease  of  Cq  with  downwind  distance  x  is 


dCp 

dx 


Cq(J  z 


1  dXT  y 

dx 


+ 


(Tzcrt 


1  da^ 
(jy(y"l  dx 


(9) 


Firstly,  note  that  the  total  concentration  flux  Qm  oc  CoCTyCr^  is  a  conserved  quantity  (more 
specifically,  for  a  Gaussian  mean  concentration  distribution  with  S  constant,  the  total  con¬ 
centration  flux  through  the  vertical  y-z  plane  is  Qm  =  2TrSCo<7y(^z)-  For  the  obstacle  array 
studied  here,  it  was  found  that  dcxy/dx  and  dajdx  (plume  growth  rates  in  the  lateral  and 
vertical  directions)  are  comparable  to  an  open-terrain  plume  dispersing  under  similar  at¬ 
mospheric  conditions.  Hence,  on  examination  of  the  preceding  equation,  it  is  seen  that  the 
slow  rate  of  decay  of  the  centreline  mean  concentration  in  the  array  plume  at  the  greater 
downwind  distances  from  the  source  arises  from  the  increased  spreads  of  the  array  plume 
in  both  the  lateral  and  vertical  directions  (due  to  the  rapid  initial  plume  growth  in  the 
“inflationary”  phase  of  plume  development). 


Concentration  fluctuation  intensity 


The  parameter  most  frequently  used  to  characterize  concentration  fluctuations  is  the  flue- 
tuation  intensity,  i  s  aJC.  where  a;,  =  [(x-CfC  i**  concentration  standard 
deviation  (x  denotes  the  instantaneous  concentration).  In  order  to  consider  the  downwind 
development  of  the  concentration  fluctuations  in  an  array  plume  in  more  detail,  we  focus 
attention  on  an  example  of  concentration  fluctuation  profiles  in  the  lateral  and  vertical 
directions  obtained  from  the  data  in  experiment  number  6  where  the  approach  flow  was 
nearly  perpendicular  to  the  front  of  the  obstacle  array.  Figure  11(a)  displays  the  lateral 
profiles  of  i  at  various  non-dimensional  downwind  distances  xi^/H  from  the  source,  whereas 
Figure  11(b)  exhibits  a  vertical  profile  of  i  along  the  mean  plume  centreline  for  this  exper¬ 
iment.  In  Figure  11(a).  the  labels  1,  2,  3,  and  4  refer  to  the  sequence  of  four  horizontal 
sampling  lines  in  the  experiment  (Figure  1)  with  lines  1  and  4  located  at  the  most  southerly 
and  most  northerly  positions,  respectively,  in  the  array. 

The  shapes  of  these  profiles  are  similar  to  those  obtained  from  plumes  dispersing  over  open 
terrain  (Mylne  and  Mason  [21];  Mylne  [22];  Yee  et  al.  [24];  Yee  et  al.  [25]).  In  particular. 
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the  lateral  cross-sections  of  fluctuation  intensity  are  seen  to  increase  towards  the  edges  of 
the  plume  from  a  minimum  at  the  plume  centreline  and  exhibit  an  approximate  bilateral 
symmetry  about  the  mean  plume  centreline  (implying  that  the  mean  wind  direction  was 
approximately  perpendicular  to  the  horizontal  sampling  lines  of  detectors).  The  vertical 
profile  varies  only  slightly  between  z/H  =  Q  and  2  (  within  the  roughness  sublayer  of  the 
obstacle  array),  but  increases  rapidly  higher  up  as  the  extreme  upper  edge  of  the  plume  is 
approached.  Despite  similarities  in  the  general  shapes  of  the  lateral  and  vertical  profiles 
of  fluctuation  intensity  to  those  observed  in  open-terrain  plumes,  the  curves  of  i  here  are 
nevertheless  wider  and  flatter,  another  manifestation  of  the  observed  increase  in  the  lateral 
and  vertical  mean  plume  spreads  (Figure  9). 

For  the  case  of  an  open-terrain  plume,  C  (the  variation  about  the  mean  concentration 
is  at  least  as  large  as  the  mean  concentration  itself)  along  the  mean-plume  centreline  {io  > 
1).  However,  from  Figure  11.  it  is  seen  that  centreline  values  of  i  are  considerably  less 
than  unity  for  plumes  dispersing  through  an  obstacle  array.  Hence,  it  appears  that  the 
decrease  in  the  concentration  standard  deviation  is  much  greater  than  even  the  reduction 
in  the  mean  concentration  in  the  array  (Figure  10),  implying  a  dramatic  reduction  in  the 
strength  of  concentration  fluctuations.  Clearly,  the  change  from  an  open-terrain  plume  to 
a  plume  dispersing  in  an  obstacle  array  has  a  very  large  effect  on  i,  reducing  it  by  a  factor 
of  between  2  and  5  along  the  plume  centreline,  depending  on  xi/H.  Quantitatively  similar 
results  were  obtained  by  Davidson  et  al.  [9]  for  the  centreline  behaviour  of  i,  where  fo  =  1.6 
was  reported  for  the  control  (open-terrain)  plume  in  contrast  to  the  significantly  smaller 
value  of  io  =  0.4  for  the  array  plume  in  an  equivalent  position. 

A  first  qualitative  explanation  for  the  drastic  reduction  of  i  in  the  array  plume  is  that  the 
lateral  and  vertical  meandering  of  the  instantaneous  plume  in  the  obstacle  array  is  reduced 
significantly  in  comparison  to  the  open-terrain  plume.  The  reduction  in  plume  meander 
within  the  array  can  probably  be  attributed  to  two  physical  mechanisms;  (1)  the  possible 
reduction  in  the  scale  of  turbulence  A  between  the  obstacles  to  that  of  the  obstacles  (~  H); 
and,  (2)  the  observed  increase  in  the  lateral  and  vertical  spreads  of  the  plume  (Figure  9). 
Hejice,  the  dramatic  increase  in  the  ratio  of  the  plume  spread  to  the  integral  turbulence  scale 
within  the  array  relative  to  its  value  in  open  terrain  would  explain  the  reduced  meandering 
of  the  array  plume,  and  therefore  a  decrease  in  the  value  of  i. 

However,  we  note  that  the  dramatic  reduction  in  i  for  plumes  within  the  array  cannot  arise 
solely  from  the  reduction  in  the  meandering  contributions.  Figure  12  compares  the  vertical 
profile  of  i  measured  within  the  obstacle  array  (Run  6)  with  a  similar  vertical  profile  of  i 
measured  at  approximately  the  same  downwind  distance  from  the  source  in  open  terrain. 
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The  data  for  the  open-terrain  plume  were  obtained  from  a  previous  experiment  (Yee  et 
al.  [25]).  In  addition,  the  vertical  profile  of  the  conditional  fluctuation  intensity  ip  for  the 
open-terrain  plume  is  also  shown  in  Figure  12,  with  the  subscript  p  used  heie  to  denote 
a  conditional  statistic  which  is  determined  from  the  non-zero  concentrations  only.  Firstly, 
because  the  intermittency  for  the  plume  within  the  array  is  small  (the  concentration  time 
series  within  the  array  plume  appears  as  a  continuous  signal  with  practically  no  periods 
where  the  concentration  was  zero)  i  ^  ip  within  the  mean-plume  core  of  the  array  plume. 
Note  that  i  for  the  array  plume  is  less  than  ip  for  the  open-terrain  plume  (Figure  12) 
throughout  most  of  the  vertical  extent  of  the  plume.  This  implies  that  the  significant 
reduction  in  i  for  the  array  plume  is  due  not  only  to  the  reduction  in  the  plume  meander, 
but  also  to  vigorous  mixing  of  material  internal  to  the  plume,  smoothing  out  the  in-plume 
concentration  fluctuations  and  transferring  the  concentration  variance  (energy)  to  smaller 
scales  through  the  inertial-convective  cascade.  In  summary,  the  decrease  in  the  plume 
meander  and  increase  in  the  dissipation  of  the  concentration  fluctuations  for  the  array 
plume  leads  to  the  observed  dramatic  decrease  in  i. 


The  observed  smoothing  of  the  internal  concentration  fluctuations  in  the  array  plume  can 
be  understood  as  follows.  Vorticity  is  generated  at  the  various  solid  surfaces  of  the  obstacles 
in  the  array  through  the  no-slip  condition.  This  vorticity  is  then  transported  throughout 
the  fluid  between  the  obstacles  by  convection  and  diffusion  and  also  will  be  enhanced  in  the 
canopy  flows  via  the  mechanism  of  vortex-stretching.  This  results  in  the  small-scale,  high- 
intensity  turbulence  that  is  characteristic  of  the  flow  within  the  obstacle  array.  The  fine-scale 
structure  of  the  canopy  turbulence  has  a  significant  effect  on  the  instantaneous  structure 
of  the  array  plume.  The  concentration  fluctuations  in  the  array  plume  are  smoothed  by 
the  vigorous  stirring  and  mixing  of  the  material  by  the  fine-scale,  high-intensity  turbulence 
in  the  canopy,  leading  to  a  drastic  reduction  in  the  length  scales  as  material  volumes  are 
thinned.  In  particular,  the  fine-scale  random  vortical  structures  in  canopy  turbulence  of 
length  scale  I  ~  a,  (o-j  is  the  mean  instantaneous  plume  width)  leads  to  the  entrainment 
or  engulfment  of  clean  air  parcels  into  the  plume,  followed  by  the  stretching  and  folding 
of  material  interfaces  (stirring)  which  generates  internal  concentration  variance  within  the 
instantaneous  plume.  This  straining  of  material  volumes  continues  until  length  scales  are 
reduced  to  where  molecular  diffusion  becomes  important.  In  the  array  plumes,  this  pro¬ 
cess  of  transferring  the  variance  to  smaller  scales  (inertial-convective  subrange  of  plume 
motions)  appears  to  be  significantly  accelerated,  rapidly  smoothing  out  the  concentration 
fluctuations. 

Figure  13  displays  the  decay  of  the  plume  centreline  fluctuation  intensity  io  with  normalized 
downstream  distance  xl/H  from  the  source  for  an  open-terrain  plume  and  for  array  plumes 
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with  various  centreline  directions  |0|  of  dispersion  through  the  array.  Note  that  the  array 
plumes  exhibit  significantly  lower  centreline  fluctuation  intensities  for  all  plume  centreline 
directions  of  dispersion  in  comparison  with  the  open-terrain  plume  at  any  given  xl/H. 
However,  the  rate  of  decrease  of  zq  in  the  array  plumes  for  xl/H  >  \s  comparable 
to  the  open-terrain  plume.  This  rate  of  decrease  is  rather  slow,  and  it  is  possible  that 
zo  asymptotically  approaches  a  non-zero  constant  value,  rather  than  decreases  to  zero.  If 
the  former  case  is  correct,  it  appears  that  zo  in  the  open-terrain  plume  reaches  a  near 
constant  value  of  about  1.0  at  long  range,  in  contrast  to  the  array  plume  where  zq  appears 
to  asymptotically  approach  a  near  constant  value  of  about  0.2  independent  of  the  direction 
of  plume  centreline  dispersion. 

The  downwind  decay  of  centreline  fluctuation  intensity  zq  stratified  by  source  height  for  a 
fixed  plume  centreline  direction  6  of  dispersion  (approximately  or  better)  and  by  the  mean 
plume  centreline  direction  of  dispersion  are  shown  in  Figures  14(a)  and  14(b),  respectively. 
We  note  that  for  sources  above  the  “urban”  canopy,  the  initial  fluctuation  intensity  is  larger 
than  for  sources  inside  the  canopy.  This  is  the  result  of  the  much  stronger  mixing  processes 
within  the  canopy  arising  from  the  small-scale,  high-intensity  turbulence  generated  by  the 
obstacles  of  the  array.  However,  the  initial  decay  rate  of  the  plume  fluctuation  intensity 
for  the  source  above  the  canopy  top  is  greater  than  that  for  the  source  within  the  canopy, 
with  the  result  that  both  the  magnitude  and  rate  of  decay  for  both  array  plumes  become 
comparable  at  greater  downwind  distances  from  the  source  (the  initial  source  conditions  are 
rapidly  “forgotten”  as  the  plume  material  is  subject  to  the  continued  stirring  and  mixing 
processes  within  the  canopy). 

There  appears  to  be  a  weak  dependence  of  the  fluctuation  intensity  level  on  the  direction 
of  plume  centreline  dispersion  through  the  array,  which  in  turn  is  a  manifestation  of  the 
mean  wind  direction  of  the  upstream  flow  (Figure  8).  In  particular,  fluctuation  intensity 
levels  in  the  array  plumes  are  larger  when  the  mean  wind  direction  is  nearly  perpendicular 
to  the  front  of  the  array  (|q|  ^  0°)  than  when  the  mean  wind  direction  is  at  some  oblique 
angle.  This  may  be  due  to  the  fact  that  when  the  wind  direction  is  nearly  perpendicular  to 
the  front  face  of  the  array,  the  obstacles  are  aligned  with  the  wind;  whereas,  when  the  wind 
direction  is  at  some  oblique  angle  to  the  front  face  of  the  array,  there  are  no  longer  channels 
through  which  the  flow  can  pass  and  the  array  behaves  more  like  a  staggered  array  than  an 
aligned  array.  The  turbulent  mixing  of  plume  material  appears  to  be  more  vigorous  in  the 
staggered  array  where  the  flow  can  impact  onto  neighboring  obstacles  than  in  the  aligned 
array  where  the  flow  is  aligned  with  the  street  canyons. 
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Concentration  probability  density  function 


The  one-point  probability  density  function  (pdf)  of  concentration  provides  information  on 
the  distribution  of  instantaneous  concentration  measured  at  a  fixed  point  in  the  plume. 
In  particular,  the  concentration  pdf  summarizes  all  the  one-point  statistics  of  the  plume 
fluctuations  (e.g.,  all  the  concentration  moments,  the  intermittency  factor,  peak-to-mean 
concentration  ratio,  etc.)  and  thereby  clearly  provides  more  complete  statistical  information 
than  any  of  the  individual  fluctuation  statistics  presented  above  (e.g.,  mean  concentration, 
fluctuation  intensity) . 

Histograms  of  the  instantaneous  concentration  x  normalized  by  the  mean  concentration  C 
were  compiled  by  sorting  the  concentration  data  into  bins  and  appropriately  normalizing  the 
bin  counts  to  form  the  concentration  pdf  f{x/C).  Figures  15  and  16  display  the  measured 
pdfs  of  x/C  along  the  mean  plume  centreline  at  various  non-dimensional  downwind  distances 
xj^jH  from  the  source  and  at  various  positions  through  a  vertical  plume  cross-section  at 
a  fixed  downwind  distance,  respectively.  Experiment  13  was  the  source  for  these  examples 
which  were  selected  in  order  to  illustrate  the  change  in  the  shape  of  the  pdf  with  downwind 
distance  from  the  source  and  with  height  above  the  ground. 

From  Figure  15,  it  is  seen  that  the  concentration  pdfs  at  the  mean  plume  centreline  at 
downwind  distances  of  xi/H  —  19.8  and  30.7  (closer  to  the  source)  exhibit  a  bimodal 
form — broad  peaks  are  observed  at  x/C"  ~  0-3  and  1.2  (arguably,  this  pdf  form  might  even  be 
interpreted  as  being  trimodal  with  the  third  broad  peak  at  x!^  ~  1.7).  Farther  downwind 
from  the  source,  these  two  peaks  in  the  concentration  pdf  merge  and  the  pdf  assumes  a 
unimodal  form  that  becomes  increasingly  concentrated  at  x/C"  ~  1-0  as  the  intermittency 
in  the  array  plume  tends  to  unity  and  the  contaminant  parcels  in  the  plume  become  more 
and  more  well  mixed.  The  lower  and  upper  tails  of  the  concentration  pdf  become  shorter 
and  shorter  with  increasing  downwind  distance  from  the  source  as  the  continued  small-scale 
mixing  in  the  array  results  in  the  homogenization  of  the  in-plume  structure. 

Figure  16  documents  the  evolution  of  the  concentration  pdf  with  height  zj H  above  the 
ground  surface.  Within  the  canopy,  the  pdf  assumes  a  unimodal  Gaussian-like  form  that  is 
concentrated  at  x/C  ~  1.0  implying  a  “well-mixed”  plume  structure  that  results  from  the 
vigorous  small-scale  stirring  and  mixing  of  the  canopy  turbulence.  Molecular  diffusion  en¬ 
hanced  by  concentration-gradient  straining  (or,  stretching  and  folding  of  material  filaments) 
arising  from  the  small-scale,  high-intensity  canopy  turbulence  produces  a  systematic  reduc¬ 
tion  of  the  central  concentration  moments  (e.g.,  mean  concentration,  concentration  variance) 
and  of  the  range  of  concentration  values  over  which  the  pdf  spans.  Above  the  canopy,  at 
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heights  between  about  zjH  =  2.0  and  3.5,  the  concentration  pdf  transitions  to  a  bimodal 
form  with  a  narrow  sharp  peak  at  x/C*  ~  1-3,  and  a  broader  weaker  peak  at  a  lower  non¬ 
zero  concentration  of  x/C"  ~  0.3.  At  still  greater  heights  above  the  canopy  (for  z/H  4.0), 
the  concentration  pdf  evolves  from  a  bimodal  form  towards  a  unimodal  form  again,  but 
now  the  unimodal  form  is  exponential-like  rather  than  Gaussian-like.  In  particular,  the 
amplitude  of  the  peak  at  x/C  >  1.0  decreases,  while  that  at  or  near  zero  concentration 
increases.  Because  the  intermittency  effect  is  reflected  in  the  probability  of  occurrence  of 
zero  or  near  zero  concentrations,  the  increase  in  the  zero  or  near  zero  concentration  peak 
in  the  pdf  simply  implies  that  there  is  a  higher  probability  of  encountering  clean  ambient 
air  parcels  near  the  upper  edge  of  the  array  plume  (and,  hence,  a  higher  probability  of 
observing  concentrations  close  to  zero). 

Interestingly,  a  comparison  of  Figure  15  with  Figure  16  shows  that  there  is  a  remarkable 
resemblance  between  the  evolution  of  the  concentration  pdf  with  increasing  height  above  the 
ground  (at  a  fixed  downwind  distance  from  the  source)  and  that  with  decreasing  downwind 
distance  from  the  source  (at  the  fixed  height  above  the  ground).  The  current  experiments  of 
dispersion  in  the  obstacle  array  did  not  measure  the  concentration  pdf  in  the  very  near-field 
region  of  the  source  (at  a  range  of  downwind  distances  within  the  first  row  or  two  of  the 
obstacle  array).  It  is  interesting  to  consider  whether  the  variation  in  the  concentration  pdf 
at  shorter  range  would  evolve  into  a  unimodal  exponential-like  form  in  analogy  with  the 
evolution  of  the  concentration  pdf  form  at  greater  heights. 

The  transition  in  the  concentration  pdf,  from  the  exponential-like  unimodal  form  higher  up 
in  the  plume  above  the  array  {z/H  >  4),  to  a  bimodal  form  with  broad  peak  concentrations 
at  x/C*  ~  0  and  x/C"  ~  1  in  the  height  range  between  about  zjH  =  2  and  3.5,  and 
finally  to  a  Gaussian-like  unimodal  form  with  the  peak  at  x/C  ~  1  for  z/H  <  1,  implies 
a  large  change  in  the  distribution  of  both  low-amplitude  events  and  high-amplitude  bursts 
of  plume  concentration  with  height.  In  particular,  the  concentration  pdfs  measured  within 
the  obstacle  array  show  shorter  lower  and  upper  tails  than  those  measured  further  up 
in  the  plume  above  the  array,  implying  that  both  very  low  and  very  high  instantaneous 
concentrations  (relative  to  C)  are  observed  with  a  smaller  frequency  within  the  canopy. 
These  features  are  consistent  with  the  general  appearance  of  the  instantaneous  concentration 
time  series  measured  simultaneously  in  a  vertical  cross-section  through  the  plume. 

Figure  17  displays  the  ratio  of  the  peak  value  of  concentration  X99  to  the  concentration 
standard  deviation  in  a  vertical  cross-section  at  the  mean  plume  centreline  for  experiment 
number  13.  Here,  the  peak  concentration  X99  is  defined  as  the  value  which  is  exceeded  only 
1%  of  the  time.  Interestingly,  approximately  constant  over  the  plume  cross-section 
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for  heights  above  the  roughness  sublayer  {zjE  >  2),  but  increases  within  the  roughness 
sublayer  above  the  canopy  and  into  the  canopy  layer.  In  particular,  above  the  roughness 
sublayer,  the  quantity  assumes  a  nearly  constant  value  of  about  4.5,  which  is  equal  to 

the  value  obtained  through  the  cross-section  of  an  open-terrain  plume  (Fackrell  and  Robins 
[15]).  However,  within  the  canopy  layer,  Xga/o’x  increases  to  a  value  between  about  6  and  7, 
implying  that  the  reduction  in  the  peak  value  of  concentration  is  less  than  the  concentration 
standard  deviation  here.  This  is  manifested  also  in  the  shape  of  the  concentration  pdf 
although  the  concentration  pdfs  measured  within  the  canopy  are  Gaussian-hke  (Figure  16), 
an  interesting  characteristic  of  these  distributions  is  their  slight  asymmetry  in  the  sense  that 
the  upper  tail  of  the  distribution  appears  to  be  slightly  elongated  compared  to  the  lower 
tail.  This  is  probably  the  result  of  the  rapid  turbulent  mixing  of  the  plume  material  by  the 
small-scale,  high-intensity  turbulence  within  the  canopy  which  effectively  removes  whatever 
large-scale  concentration  gradients  there  are  in  the  array  plume,  forming  patches  of  near 
constant  concentration  which  can  conform  with  the  existing  large-scale  gradients  only  by 
the  formation  of  “fronts”  (sheets  of  intense  concentration  gradient).  This  interpretation 
appears  consistent  with  the  appearance  of  the  concentration  time  traces  measured  within 
the  array,  where  "plateau  and  cliff”  structures  are  evident. 


It  is  useful  to  develop  a  mathematical  model  for  the  concentration  pdf  of  the  array  plume. 
To  this  end,  it  is  pragmatic  to  investigate  pdf  forms  which  can  be  fully  specified,  including 
the  intermittency  effect,  using  only  two  parameters.  Three  commonly  used  pdf  models  for 
the  concentration  which  satisfy  this  constraint  are  the  exponential  pdf  (Barry  [36]),  the 
clipped  normal  pdf  (Lewellen  and  Sykes  [37]),  and  the  clipped  gamma  pdf  (Yee  and  Chan 
[38]).  The  exponential  pdf  assumes  the  form 

/(x)  =  ^exp  (^-^)  +  (1 -7)<^(x)>  (10) 

where  7  is  the  intermittency  factor,  defined  as  the  probability  that  the  concentration  is 
non-zero  and  6{-)  is  the  Dirac  delta  function.  The  clipped  normal  pdf  is  given  by 


where 


7  =  7(7‘: 


l^c 

v/2(Tc 


Here,  /ic  and  CTc  are  location  and  scale  parameters,  respectively  and  erf(-)  is  the  error 
function.  Unlike  a  Gaussian  pdf,  f^c  C  for  the  clipped  normal  pdf.  Finally,  the  clipped 


gamma  pdf  is  defined  as 


exp(-(x  +  A)/s) 
sT{k) 


+  (1  -7)(5(x), 


(13) 
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where  s,  k,  and  A  are  scale,  shape,  and  shift  parameters,  respectively,  that  are  constrained 
as  follows: 

7  =  =  (14) 


where  {x/CY  is  the  normalized  mean-square  concentration.  Here,  r(x)and  r(jz;a;)  denote 
the  Gamma  function  and  incomplete  Gamma  function,  respectively.  We  note  that  for 
the  case  where  7  — >  1“  [or,  equivalently  i  =  {{x/CY  +  l)  iii  accordance  to 

Equation  (14)],  A  — >  O"*",  and  the  clipped  gamma  distribution  reduces  to  the  standard 
gamma  distribution.  Because  7  =  1  is  the  upper  bound  for  the  intermittency  factor,  for 
cases  where  i  <  the  standard  gamma  distribution  with  parameters  s  and  k  is  applied 
(A  =  0  here). 

To  determine  how  well  the  model  pdfs  represent  the  concentration  data,  it  is  useful  to 
construct  quantile-quantile  (Q-Q)  plots  that  compare  the  order  statistics  from  the  sample 
data  with  those  from  the  model  distributions.  Recall  that  the  quantile  of  order  g  (0  < 
g  <  1),  designated  by  Xioogi  is  the  unique  solution  of  the  equation  ^(xioo?)  =  Pr(x  < 
Xioog)  =  Q  (the  concentration  level,  Xioogi  exceeds  the  instantaneous  concentration  values 
with  probability  q). 

Figure  18  presents  the  Q-Q  plots  of  the  observed  data  quantiles  against  the  predicted 
(model)  quantiles  for  2  of  the  3  theoretical  distributions  summarized  above.  We  note  that 
the  exponential  distribution  cannot  be  used  to  model  the  plume  concentrations  in  the  exam¬ 
ples  of  Figure  18  because  the  observed  fluctuation  intensities  were  less  than  unity.  The  data 
are  taken  at  various  heights  along  the  centreline  of  an  array  plume  at  a  fixed  downwind  dis¬ 
tance  from  the  source  (Trial  13).  An  exact  fit  between  the  observed  and  model  distribution 
results  in  a  straight  line  with  slope  1  and  intercept  0  when  the  abscissa  (observed  quantiles) 
and  ordinate  (model  quantiles)  have  equal  scales.  In  the  interpretation  of  the  Q-Q  plots, 
one  should  keep  in  mind  that  departures  of  the  measured  quantiles  from  the  straight  line 
in  the  extreme  lower  and  upper  tails  of  the  distribution  may  simply  reflect  a  lack  of  statis¬ 
tical  convergence  (viz.,  inadequate  statistical  sampling)  in  these  extreme  tails.  With  this 
caveat,  it  appears  that  the  clipped  normal  and  clipped  gamma  distributions  can  represent 
adequately  the  observed  concentration  pdf  over  much  of  its  entire  range.  Nevertheless,  the 
clipped  gamma  distribution  appears  to  be  a  better  representation  of  the  observed  distribu¬ 
tion  than  the  clipped  normal  distribution.  In  particular,  both  the  upper  and  lower  tails  of 
the  observed  distribution  are  better  characterized  by  the  clipped  gamma  distribution.  The 
probability  of  non-zero  concentrations,  7,  is  underpredicted  by  the  clipped  normal  distribu¬ 
tion  with  the  consequence  that  the  lower  tail  of  the  clipped  normal  pdf  is  generally  lighter 
(shorter)  than  that  for  the  observed  concentration  pdf.  This  tendency  of  the  clipped  normal 
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pdf  to  underpredict  the  lower  tail  of  the  concentration  distribution  has  been  noted  earlier 
by  Yee  and  Chan  [38]  for  open-terrain  plumes.  Finally,  the  clipped  gamma  pdf  appears 
to  predict  the  slight  asymmetr}'  in  the  measured  in-array  concentration  pdfs  as  well  as  the 
elongation  of  the  observed  upper  tail  better  than  the  clipped  normal  pdf. 


Concentration  power  spectra 


The  normalized  power  spectral  density  of  concentration  fluctuations,  nSy{n)/o'^,  plotted 
against  frequency,  n,  are  exhibited  in  Figure  19  for  a  number  of  different  heights  above 
the  ground  at  the  centreline  of  an  array  plume  at  a  range  of  around  xl/H  =  30.7  (ex¬ 
periment  13).  Here,  5x(n)  is  the  concentration  power  spectral  density.  The  concentration 
power  spectral  density  was  estimated  by  dividing  the  entire  time  series  into  a  number  of 
blocks,  applying  a  Parzen  window  or  taper  to  each  block  to  suppress  end  effects  due  to 
finite  truncation,  computing  the  Fast  Fourier  Transform  (FFT)  of  each  block,  and  then 
averaging  the  squared  modulus  of  the  FFT  (the  raw  periodogram)  at  each  frequency  over 
the  different  blocks  to  obtain  an  ensemble  average  estimate  of  the  power  spectral  density 
at  this  frequency.  The  spectra  in  Figure  19  were  measured  on  the  32-m  tower  near  the 
centre  of  the  obstacle  array.  The  area  under  the  various  spectra  in  Figure  19  have  been 
equalized  by  normalization  by  the  variance  of  the  concentration  fiuctuations  in  accordance 
to  the  relationship 

rOO 

/  Sy.{n)dn  —  (15) 

Jo 

Hence,  the  concentration  power  spectra  shown  in  Figure  19  are  one-sided  spectra. 

For  concentration  spectra  measured  above  the  canopy  top,  there  appears  to  be  very  little 
variation  in  the  form  of  the  spectra.  All  these  spectra  exhibit  a  pronounced  peak  at  about  0.4 
Hz,  marking  the  beginning  of  an  inertial-convective  subrange.  The  establishment  of  small- 
scale  motion  in  the  plume  concentration  fluctuations  is  often  connected  with  the  existence  of 
an  inertial-convective  subrange  in  which  the  slope  of  log(S';^(n))  versus  log(n)  is  —5/3  (Mylne 
and  Mason  [21]).  The  range  of  scales  associated  with  the  inertial-convective  subrange 
corresponds  to  eddies  on  smaller  scales  than  the  instantaneous  plume,  and  these  eddy  scales 
are  responsible  primarily  for  the  internal  mixing  process  which  transfers  fluctuation  variance 
from  large  scale  to  the  very  small  scales  where  it  can  be  dissipated  by  molecular  diffusion.  It 
is  seen  that  concentration  spectra  measured  above  the  canopy  exhibit  a  moderate  frequency 
range  in  which  the  slope  is  of  the  order  -5/3  (or,  equivalently,  nS^{n)  ~  n-^/^)  with 
the  upper  cutoff  being  at  about  4  Hz.  Hence,  these  spectra  exhibit  an  inertial-convective 
subrange  spanning  about  one  decade  in  frequency  between  about  0.4  Hz  and  4  Hz  associated 
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with  the  range  of  scales  at  which  concentration  variance  (energy)  is  transferred  from  large 
to  small  scales  with  negligible  molecular  dissipation  and / or  interaction  with  the  gradients  in 
the  mean  concentration  field  (i.e.,  gradient  production).  This  range  is  associated  with  those 
frequencies  that  correspond  to  scales  of  motion  that  are  small  relative  to  the  instantaneous 
plume  width. 

In  contrast,  the  power  spectra  measured  in  the  canopy  layer  {zlH  <  1)  are  significantly 
different  than  those  measured  above  the  canopy  top.  In  particular,  a  rapid  drop-off  in  the 
concentration  power  spectra  occurs  in  a  narrow  band  of  frequencies  just  greater  than  the 
spectral  peak  at  about  0.4  Hz  (i.e.,  0.4  Hz  <  n  <  4  Hz)  where  the  spectral  slope  is  seen  to 
be  significantly  more  negative  than  -5/3.  This  drop-off  in  the  spectral  density  is  greater 
for  positions  that  are  deeper  in  the  canopy  (closer  to  the  ground  surface).  In  other  words, 
there  appears  to  be  a  monotonic  reduction  in  spectral  level  for  n  >  0.4  Hz  with  decreasing 
height  above  the  ground  within  the  urban  canopy  layer.  At  frequencies  above  about  4  Hz, 
all  spectra  (measured  both  within  and  above  the  urban  canopy)  are  seen  to  drop-off  rapidly 
with  a  logarithmic  slope  that  is  significantly  more  negative  than  —5/3. 

The  significant  “loss”  of  spectral  power  at  frequencies  n  slightly  greater  than  the  spectral 
peak  at  about  0.4  Hz  for  in-canopy  concentration  power  spectra  can  be  attributed  to  the 
enhanced  small-scale  stirring  and  mixing  of  plume  material  within  the  canopy.  Indeed, 
the  presence  of  the  solid  surfaces  of  the  obstacles  in  the  array  leads  to  the  generation  of 
fine-scale  turbulence  or  small-scale  random  vortical  structures  in  the  wake  flows  that  act 
to  mix  the  array  plume  elements  internally,  smoothing  out  the  concentration  fluctuations. 
This  action  is  expected  to  result  in  a  greater  concentration  variance  loss  from  the  high- 
frequency  (small-scale)  end,  with  the  result  that  a  greater  fraction  of  the  total  variance  of 
in-canopy  concentration  spectra  is  concentrated  at  the  low  frequency  end  compared  with 
concentration  spectra  measured  above  the  canopy.  Hence,  it  appears  that  the  fine-scale 
wake  turbulence  results  in  a  rapid  and  efficient  mixing  of  the  plume  material  within  the 
canopy  layer,  resulting  in  a  short  circuit  of  the  eddy  cascade  of  concentration  variance 
(energy)  through  the  inertial-convective  subrange  of  scales. 

Interestingly,  the  location  of  the  spectral  peak  in  the  plots  of  nS^{n)  versus  log(n)  does  not 
appear  to  change  as  we  descend  into  the  “urban”  canopy.  This  spectral  peak  in  the  array 
plume  concentration  fluctuations  appears  to  remain  invariant  with  increasing  downwind 
distance  from  the  source.  In  particular,  the  location  of  the  concentration  spectral  peak  Up 
is  at  about  0.4  Hz;  or,  equivalently,  UpH/Sn  ^  0.51  ±  0.02  (where  Sh  denotes  the  mean 
horizontal  wind  speed  measured  at  the  canopy  top  zjH  =  1).  This  behaviour  contrasts 
with  the  case  of  an  open-terrain  plume  where  it  was  generally  found  that  the  spectral 
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peak  frequency  shifts  toward  lower  frequencies  with  increasing  downwind  distance  from  the 
source  (Mylne  and  Mason  [21];  Yee  et  al.  [23];  Yee  et  al.  [24]).  In  particular,  for  an 
open-terrain  plume  Up  is  proportional  to  S/oy  (Mylne  and  Mason  [21]),  where  S  is  the 
mean  horizontal  wind  speed  and  ay  is  the  lateral  plume  spread,  whereas  in  the  array  plume 
UpH/Sn  ~  0-51  •  Hence,  the  spectral  peak  in  the  array  plume  concentration  fluctuations 
appears  to  be  imposed  by  the  topological  characteristics  of  the  obstacle  array,  and  appears 
to  be  related  to  the  time  scale  H/Sh  of  turbulent  eddies  of  the  scale  of  the  obstacles  (wake 
scales).  Furthermore,  the  time  scale  H/Sh  is  also  comparable  to  the  time  it  takes  plume 
material  entrained  into  the  wake  of  an  obstacle  of  height  H  to  diffuse  out  of  it  (Vincent 
[39]).  In  summary  then,  the  frequency  of  the  spectral  peak  which  does  not  change  as  we 
descend  into  the  canopy  can  be  interpreted  as  an  inverse  time  scale  for  the  energy-containing 
concentration  eddies  in  the  array  plume,  eddies  that  appear  to  be  related  to  the  wake  scales 

of  motion  in  the  canopy. 

Figure  20  shows  an  example  of  the  cross-spectral  density  of  array  plume  concentration  fluc¬ 
tuations  (which  has  been  decomposed  into  the  squared  coherency  nl^in)  and  phase  4>i2{n)) 
between  two  points  in  the  plume  measured  below  the  canopy  height.  The  squared  coherency 
Ki2(n)  between  concentration  time  series  xi(^)  and  X2{t)  measured  at  z  =  zi  and  z  =  Z2, 
respectively,  at  frequency  n  is  the  squared  modulus  of  the  cross-spectral  density  S'i2(n) 
at  n,  normalized  by  the  products  of  the  two  auto-spectral  densities  5i(n)  and  S2{n),  so 
Ki2(n)  =  |5'l2(n)|V(‘5l(«)‘S'2(^^))•  The  phase  (I>i2{n)  of  xi(0  and  X2{t)  at  frequency  n  is 
the  argument  (angle)  of  the  cross-spectra  density  Suin).  Together,  «i2(n)  and  (pnin)  de¬ 
scribe  in  the  frequency  domain  the  correlation  structure  between  the  concentration  processes 
measured  at  zjH  =  0.4  (point  1)  and  zjH  =  0.8  (point  2).  The  phase  spectrum  0i2(n) 
represents  the  “average  value”  of  the  phase  shift  between  the  components  in  xi^/H  =  0.4) 
and  x(^/-H  =  0.8)  at  frequency  n  (defined  mod  27r)  and  is  conventionally  taken  to  be  in 
the  interval  (-tt.tt).  However,  the  phase  spectrum  defined  on  the  interval  (-7r,7r)  exhibits 
discontinuities,  and  in  order  to  avoid  this  complication,  the  phase  spectrum  displayed  in 
Figure  20  has  been  phase  “unwrapped”  to  provide  a  continuous  function  of  n. 

The  squared  coherency  and  phase  spectra  exhibited  in  Figure  20  show  a  sudden  change 
from  high-coherence,  in-phase  variations  for  n  <  0.4  Hz  to  incoherent  random  fluctuations 
for  n  >  0.4  Hz.  The  correlation  coefficient  between  the  two  concentration  time  series 
was  0.79  for  this  example  (taken  from  Run  13)  and  can  be  seen  to  arise  almost  entirely 
from  frequencies  less  than  about  0.4  Hz  (which  coincides  also  with  the  concentration  power 
spectral  peak).  The  coherency  in  the  plume  concentration  fluctuations  measured  at  different 
heights  in  the  canopy  appear  to  be  related  to  the  wake  scales  of  motion  which  extend  to 
the  top  of  the  canopy  (and  are  determined  as  such  by  the  height  H  of  the  obstacles). 
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time  scale  variations  must  be  the  result  of  spatial  variations  in  the  plume  itself.  In  particular, 
the  increased  time  and  length  concentration  scales  within  the  canopy  are  consistent  with  the 
observation  that  the  plume  structure  is  smoother  and  less  fragmented  within  the  canopy  as  a 
consequence  of  the  vigorous  small-scale  mixing  imposed  by  the  canopy  turbulence.  Finally, 
the  integral  and  Taylor  length  scales  at  mid-canopy  height  are  comparable  to  characteristic 
dimensions  of  the  individual  canopy  elements — A^{z/H  =  0.5)  si  Lp  and  At{z/H  =  0.5)  ^ 
H,  where  Lp  =  VWH  is  a  characteristic  dimension  of  the  front  face  of  the  obstacle. 

The  Taylor  time  scale  appears  to  be  roughly  constant  over  much  of  the  vertical  plume 
cross-section,  with  a  moderate  increase  seen  to  occur  only  below  the  top  of  the  canopy.  The 
variations  in  the  Taylor  length  scale  seen  in  the  height  range  z/H  >  2  are  due  solely  to 
variations  in  the  mean  wind  speed  profile  (i.e.,  the  increase  in  the  Taylor  length  scale  with 
height  above  the  canopy  merely  reflects  the  increase  in  the  mean  wind  speed  with  height 
here) . 

The  enhanced  and  rapid  mixing  of  plume  material  in  the  obstacle  array  merges  small-scale 
plume  concentration  structures,  leading  to  a  reduction  in  the  frequency  and  number  of  low- 
amplitude  events  and  high-amplitude  bursts,  which  increases  the  concentration  fluctuation 
integral  and  Taylor  microscale  time  and  length  scales  in  the  canopy.  Hence,  canopy  turbu¬ 
lence  in  the  obstacle  array  appears  to  function  as  an  efficient  “mixer”,  reducing  the  plume 
fluctuation  intensity  and  increasing  the  plume  time  and  length  scales.  In  effect,  the  current 
measurements  suggest  that  canopy  turbulence  acts  as  a  truly  effective  mixing  flow  (flow  that 
results  in  rapid  mixing  of  material  released  into  it).  The  rapid  initial  reduction  in  the  plume 
fluctuation  intensity  within  the  obstacle  array  implies  plume  material  filaments  injected  at 
any  position  in  the  “urban"  canopy  appear  to  be  stretched  (at  least  initially)  at  a  geomet¬ 
rical  or  exponential  rate,  resulting  in  greatly  enhanced  concentration  gradients  (generated 
by  the  ever  thinning  and  convoluted  filaments)  that  facilitate  rapid  mixing  by  molecular 
diffusion.  The  increased  time  and  length  scales  of  plume  concentration  fluctuations  in  the 
canopy  imply  that  the  material  distribution  in  the  instantaneous  plume  is  more  coherent 
here.  The  coherent  nature  of  the  array  plume  has  an  important  implication  for  plumes 
composed  of  inflammable  material — in  a  coherent  plume,  there  is  a  higher  probability  that 
an  ignition  source  at  a  given  position  and  time  can  result  in  complete  burn-up  because  in 
this  case,  there  is  a  greater  chance  that  there  is  a  continuous  path  of  inflammable  material 
extending  back  to  the  source. 
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Conclusions 


The  primary  purpose  of  this  paper  is  to  present  a  detailed  picture  of  the  behaviour  of 
concentration  fluctuation  statistics  in  plumes  from  point  sources  dispersing  through  a  large 
array  of  obstacles.  The  aim  is  mainly  to  provide  valuable  basic  data  on  the  structure  and 
development  of  array  plumes  that  can  be  used  to  guide  the  construction  and  assessment 
of  predictive  models  for  dispersion  in  a  built-up  area.  However,  the  basic  field  experiment 
conducted  here  has  its  own  merit  in  that  the  effect  of  a  large  group  of  obstacles  on  the 
dispersion  of  a  passive  scalar  can  be  investigated  by  analyzing  the  data  presented  herein. 
The  main  conclusions  of  this  work  may  be  summarized  as  follows. 

(1)  For  plumes  released  from  a  small  point  source  where  the  initial  plume  size  is  small 
relative  to  the  size  of  the  individual  obstacles  in  the  array,  the  lateral  and  vertical  plume 
spreads  are  increased  initially  rather  significantly  (by  a  factor  of  2  to  4)  relative  to  an 
open- terrain  plume  measured  under  similar  conditions  (“inflationary  phase”  of  array  mean 
plume  development).  However,  the  rapid  initial  increase  in  the  array  plume  dimensions  is 
followed  subsequently  by  a  rate  of  growth  that  is  approximately  comparable  to  that  of  an 
open-terrain  plume  under  the  same  atmospheric  conditions.  The  observed  initial  increase 
in  the  lateral  plume  width  was  due  to  the  large  width-to-height  {W/H  ^  4.8)  aspect  ratio 
of  the  obstacles  used  in  this  study. 

(2)  In  the  array  of  wide  and  low  obstacles  studied  here,  it  was  found  that  the  mean  plume 
centreline  concentration  in  the  array  is  generally  significantly  smaller  (by  a  factor  of  5  to 
10)  than  that  observed  in  an  open-terrain  plume  under  similar  atmospheric  conditions. 
This  observation  diff'ers  markedly  from  the  case  of  an  array  of  cubes  (unity  aspect  ratio 
obstacles)  where  it  was  found  that  the  centreline  mean  concentration  (at  the  source  height) 
was  generally  similar  to  or  greater  than  that  in  the  open-terrain  plume  (Davidson  et  al. 
[10];  Macdonald  et  al.  [11]).  Hence,  for  an  array  composed  of  obstacles  with  significant 
spanwise  aspect  ratio,  it  appears  that  the  reduction  in  mean  velocity  within  the  array 
cannot  counteract  the  significant  increases  in  both  the  lateral  and  vertical  plume  spreads. 
In  addition,  we  found  that  the  rate  of  decay  of  the  centreline  mean  concentration  in  the 
array  plume  at  downwind  distances  x^/H  >  30  was  slower  than  that  in  an  open-terrain 
plume  dispersing  under  similar  atmospheric  conditions. 

(3)  Lateral  profiles  of  the  mean  concentration  are  represented  well  with  a  Gaussian  form. 
The  crosswind  mean  concentration  profiles  in  the  array  plume  at  different  downwind  dis¬ 
tances  from  the  source  can  well  be  collapsed  onto  a  single  Gaussian  form  when  they  are 


DRDC  Sufficld  TR  2003-097 


37 


normalized  by  the  local  peak  mean  concentration  Coy  and  the  local  lateral  plume  width  ay. 
However,  vertical  profiles  of  the  mean  concentration  develop  in  a  rather  complex  manner 
with  increasing  downwind  distance  from  the  source,  and  generally  it  was  found  that  the  re¬ 
flected  Gaussian  form  is  not  always  an  ideal  description  of  the  vertical  mean  concentration 

distribution. 


(4)  The  increased  lateral  and  vertical  plume  spreads  coupled  with  the  reduction  of  the  scale 
of  turbulence  in  the  obstacle  array  results  in  a  significantly  reduced  plume  meander.  The 
reduction  in  plume  meander  and  the  smoothing  out  of  the  in-plume  concentration  fluctu¬ 
ations  results  in  a  drastic  reduction  in  the  fluctuation  intensity  of  an  array  plume  relative 
to  an  open-terrain  plume  under  similar  conditions.  In  particular,  centreline  fluctuation  in¬ 
tensities  in  the  array  plume  are  generally  a  factor  of  between  2  and  5  smaller  than  those  m 
the  open-terrain  plume  depending  on  the  downwind  distance  from  the  source.  In  addition, 
the  initial  rate  of  decay  of  concentration  fluctuations  is  increased  in  the  array  plume  in 
comparison  with  the  open-terrain  plume,  although  at  greater  downwind  distances  from  the 
source  {xl/H  >  30)  the  rate  of  decay  of  concentration  fluctuations  in  the  array  plume  ap¬ 
pears  to  be  comparable  to  that  in  the  open-terrain  plume.  The  reduction  m  the  fluctuation 
intensities  is  observed  also  in  the  array  plume  above  the  top  of  the  “urban”  canopy  and  is  a 
direct  consequence  of  the  increased  mean  wind  shear  produced  near  the  canopy  height  (e.g., 
near  zjH  ~  1,  dSjdz  1),  resulting  initially  in  the  rapid  amplification  of  scalar  gradients, 
followed  by  rapid  molecular  diffusion  across  the  steep  concentration  gradients  established 
by  the  increased  shearing  action. 


(5)  Rapid  mixing  of  plume  material  results  in  the  shortening  of  the  lower  and  upper  tails 
of  the  concentration  pdf  measured  within  the  obstacle  array,  implying  that  both  very  low 
and  very  high  concentrations  are  observed  with  smaller  frequency  here.  The  shortening  of 
the  lower  tail  in  the  concentration  pdf  is  manifested  in  the  increase  in  the  mtermittency 


factor  7  (i.e.,  7  ~  1)  for  plume  concentration  fluctuations  in  the  “urban”  canopy,  implying 
there  are  significantly  fewer  pockets  of  clean  (unmixed)  air  in  the  plume  here.  The  pdf 
of  the  array  plume  concentrations  obtained  within  the  “urban”  canopy  rapidly  develops  a 
tendency  towards  a  Gaussian-like  form,  albeit  with  a  small  asymmetry  in  the  sense  that  the 
upper  tail  of  the  distribution  is  slightly  longer  than  the  lower  tail.  The  rapid  relaxation  of 
the  concentration  pdf  to  a  Gaussian-like  form  implies  that  (where  Tc  and  are  the 

time  scales  characterizing  the  creation  and  destruction  of  plume  concentration  fluctuations). 


in  turn  implying  a  regime  where  molecular  diffusion  dominates.  In  this  regime,  a  Gaussian 
pdf  might  be  expected.  Certainly,  the  concentration  pdf  forms  observed  are  consistent  with 
the  notion  of  the  importance  of  the  continuous  operation  of  fine-scale  canopy  turbulence  m 
reducing  gradients  of  plume  concentration  and  promoting  the  blending  of  finely  laminated 
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concentration  patterns.  Finally,  we  found  that  some  simple  two-parameter  model  distribu¬ 
tions  such  as  the  clipped  normal  and  clipped  gamma  pdfs  were  capable  of  modelling  the 
distribution  of  scalar  concentrations  over  a  wide  range  of  positions  in  the  array  plume. 

(6)  The  concentration  power  spectra  measured  above  the  obstacle  array  exhibited  a  mod¬ 
erate  inertial-convective  subrange  spanning  about  a  decade  in  frequency.  There  was  a 
dramatic  reduction  in  the  spectral  power  level  in  the  inertial-convective  subrange  for  the 
part  of  the  plume  lying  within  the  obstacle  array.  Furthermore,  the  spectra  here  had  more 
of  the  total  concentration  variance  concentrated  at  the  lower  frequencies  (n  <  rip  ^  0.4 
Hz).  This  suggests  that  the  small-scale,  high-intensity  random  motions  in  the  canopy  tur¬ 
bulence  must  be  relatively  more  effective  at  smearing  out  concentration  fluctuations  at  the 
small  scales  characteristic  of  the  inertial-convective  subrange  (associated  with  the  internal 
structure  or  patchiness  in  the  plume)  than  they  are  at  the  larger  scales  characteristic  of  the 
energetic  subrange. 
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Figure  1.  A  schematic  diagram  of  the  geometry  of  the  Mock  Urban  Setting  Trial  (MUST)  obstacle  array. 
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Figure  2.  Spatial  variation  of  the  normalized  mean  horizontal  wind  speed  at 
various  heights  from  the  upwind  fetch,  through  the  obstacle  array,  and  into  the 
downwind  fetch  (a)  for  near  perpendicular  flow  and  (b)  for  a  large  obliquity  of  flow 
incidence.  The  letters  A,  B,  C,  or  D  associated  with  the  data  points  indicate  that 
the  corresponding  point  was  obtained  from  a  sonic  anemometer  on  towers  A,  B, 
C,  or  D,  respectively  (Figure  1).  The  label  “in  gap”  refers  to  a  measurement  of  the 
normalized  mean  horizontal  wind  speed  made  by  the  sonic  anemometer  (VX 
probe)  positioned  below  the  urban  canopy  height  in  the  spanwise  gap  between 
obstacles  G6  and  G7. 
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Figure  3.  Spatial  variation  of  the  normalized  turbulence  kinetic  energy  at  various 
heights  from  the  upwind  fetch,  through  the  obstacle  array,  and  into  the 
downwind  fetch  (a)  for  near  perpendicular  flow  and  (b)  for  a  large  flow 
incidence  wind  angle.  The  letters  A,  B,  C,  or  D  associated  with  the  data  points 
indicate  that  the  corresponding  point  was  obtained  from  a  sonic  anemometer  on 
towers  A,  B,  C,  or  D,  respectively  (Figure  1).  The  label  “in  gap”  refers  to  a 
measurement  of  the  normalized  mean  horizontal  wind  speed  made  by  the  sonic 
anemometer  (VX  probe)  positioned  below  the  urban  canopy  height  in  the 
spanwise  gap  between  obstacles  G6  and  G7. 


DRDC  Suffield  TR  2003-097 


47 


Figure  4.  Variation  of  the  buik  drag  coefficient  with  the  stability  parameter 
obtained  from  the  3-D  sonic  anemometer  at  the  4-m  level  on  the  32-m  vertical 
tower  near  the  centre  of  the  obstacle  array. 
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Figure  5.  Crosswind  (lateral)  mean  concentration  profiles  measured  at  two 
normalized  downwind  distances  from  the  source  within  the  obstacle  array  for 
near  perpendicular  flow  incidence.  The  data  are  extracted  from  two  horizontal 
sampling  lines  in  Trial  6.  Superimposed  on  these  profiles  are  the  best-fit 
Gaussian  distributions  (see  Eq.  (6)  for  a  definition  of  the  fitted  parameters).  The 
arrows  indicate  the  lateral  position  of  the  source. 
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Figure  6.  Vertical  profiles  of  the  mean  concentration  at  two  normalized 
downwind  distances  and  source  release  heights  for  array  plumes.  Figure  6(a) 
shows  experiment  number  6  at  a  normalized  downwind  distance  of  18.0  with  the 
source  at  a  normalized  height  of  0.06,  and  Figure  6(b)  shows  experiment 
number  10  at  a  normalized  downwind  distance  of  25.4  with  the  source  at  a 
normalized  height  of  0.7.  Superimposed  on  these  profiles  are  the  best-fit 
reflected  Gaussian  distributions  (see  Eq.  (7)  for  a  definition  of  the  fitted 
parameters). 
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Figure  7.  Comparison  of  the  vertical  normalized  mean  concentration  profiles 
along  the  plume  centreline  at  a  normalized  downwind  distance  of  48.0  from  the 
source  for  point  sources  location  1  m  (Trial  18)  and  24  m  (Trial  17)  upwind  of 
the  front  face  (atx=  0)  of  the  obstacle  array.  For  both  releases,  the  incident  wind 
direction  was  39  degrees  and  the  source  was  at  a  normalized  height  of  0.5. 
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Figure  9.  The  (a)  lateral  and  (b)  vertical  growth  of  a  plume  dispersing  in  the 
obstacle  array  with  nearly  perpendicular  flow  incidence,  compared  with  standard 
Pasquill-Gifford  dispersion  curves  for  open  terrain  for  various  atmospheric 
stability  classes  B,  C,  D,  and  E. 
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Figure  10.  The  variation  of  the  normalized  mean  plume  centreline  concentration 
with  normalized  downwind  distance  from  the  source  through  the  obstacle  array 
for  various  normalized  source  heights  and  directions  of  the  plume  centreline 
dispersion.  The  variation  is  compared  with  standard  Gaussian  plume  model 
results  using  Pasguill-Gifford  dispersion  parameters  for  open  (rural)  terrain. 
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Figure  11.  (a)  Crosswind  cross-sections  along  four  horizontal  sampling  lines 
(labeled  1,  2,  3,  and  4)  at  various  normalized  downwind  distances  from  the 
source  and  (b)  a  vertical  profile  obtained  on  the  32-m  vertical  tower  of  the 
fluctuation  intensity  in  an  array  plume  with  near  perpendicular  flow  incidence. 
Crosswind  profiles  of  fluctuation  intensity  were  measured  at  normalized 
downwind  distances  from  the  source  of  6.6,  18,  31,  and  42.  The  dashed  line 
indicates  the  lateral  position  of  the  source.  The  data  for  this  example  were 
extracted  from  Trial  6. 
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Figure  12.  Comparisons  of  vertical  profiles  of  fluctuation  intensity  measured  in  the 
array  plume  with  those  measured  in  an  open-terrain  plume  under  similar 
atmospheric  conditions  and  at  about  the  same  downwind  distance  from  the 
source.  The  vertical  profile  of  the  conditional  fluctuation  intensity  in  the  open- 
terrain  plume  is  also  shown  for  comparison. 
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Figure  13.  Plume  centreline  values  of  fluctuation  intensity  plotted  against 
normalized  downwind  distance  from  the  source  for  dispersion  over  open  terrain 
and  inside  the  obstacle  array  for  various  directions  of  plume  centreline  dispersion. 
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Figure  15.  Probability  density  function  of  normalized  concentration  measured  along 
the  centreline  of  a  dispersing  plume  within  the  array  at  a  fixed  normalized  height  of 
0.71  at  various  non-dimensional  downwind  distances  from  the  source.  The  data 
for  this  example  were  extracted  from  Trial  13. 
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Figure  16.  Vertical  profiles  of  the  probability  density  function  of  normalized 
concentration  measured  along  the  centreline  of  a  dispersing  plume  within  the 
array  at  a  fixed  non-dimensional  downwind  distance  of  30.7  from  the  source.  The 
data  for  this  example  was  taken  from  Trial  13. 
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Figure  17.  Vertical  profiles  of  the  ratio  of  the  peak  concentration  to  the 
concentration  standard  deviation  measured  along  the  centreline  of  an  array  plume 
at  a  normalized  downwind  distance  from  the  source  of  30. 7.  The  data  for  this 
example  were  extracted  from  Trial  13. 
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Figure  18.  Quantile-quantile  (Q-Q)  plots  comparing  the  normalized  concentration 
data  quantiles  with  the  associated  model  quantiles  of  the  fitted  clipped  normal 
and  clipped  gamma  distributions.  The  data  were  extracted  from  Trial  13  at 
various  heights  in  the  plume  at  a  non-dimensional  downwind  distance  from  the 
source  of  30. 7. 
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Figure  19.  Vertical  variation  of  the  normalized  power  spectrum  of  concentration 
fluctuations.  A  straight  line  with  a  slope  of -2/3  has  been  included  in  the 
logarithmic  plots  for  comparison  with  the  spectra.  The  data  were  extracted  from 
Trial  13  at  various  heights  in  the  plume  at  a  non-dimensional  downwind  distance 
from  the  source  of  30. 7. 
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Figure  20.  (a)  Squared  coherency  and  (b)  phase  spectra  between  two 
concentration  time  series  measured  at  two  different  heights  below  the  top  of  the 
obstacie  array.  The  data  for  this  example  were  extracted  from  Trial  13. 
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Figure  21.  Vertical  profiles  of  (a)  concentration  integral  time  and  length  scales 
and  (b)  Taylor  microscale  time  and  length  scales  measured  in  the  array  plume 
along  the  mean  plume  centreline  at  a  normalized  downwind  distance  from  the 
source  of  30.7.  The  data  for  this  example  were  extracted  from  Trial  13. 
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